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SYITOPBJG 


This thesis presents a study of nuclear magnetic 
relaxation in coupled spin systems in liq_uids, by steady- 
state homonuclear double resonance. The n.m.r. spectra 
of most liq_uids under high resoltition show a multiplet 
structure. These spectra normally exhibit little informa- 
tion on spin relaxation processes operative in the system. 
The relaxation in these systems can be studied by double 
resonance experiments, in which one of the transitions in 
the spectr-um is irradiated by a strong r. f. field while 
the remaining part of the spectrum is observed by a second 
weaJc r.f. field. Irradiation causes features in the 
spectrum which depend on relaxation. An analysis of 
these features loads to significant information about 
the mechanisms of relaxation. Double resonance experiments 
in this work, were performed on the coupled spin, systems 
of (i) ring protons of 2, 6-dibromoonilinc (AB 2 spin 
system) and (ii) of 2-bronotM azole (AB spin system). 

These experiments were performed with a view to study 
the systematica of the method oijd details of the 
relaxation processes. The motivation of this study and 
a brief reviov; of earlier work is given in the introductory 
chapter of the thesis. 



The notation and the density~matrix method of 
analysis, used in this work, are presented in Chapter IIA» 
The inhomogeneity of the static magnetic field is found 
to play an important role in the interpretation of the 
relaxation effects and a method of including it in the 
density-matrix calculations and of determining its 
contribution to single and double resonance transitions 
is given in Sec. IIB, The. distinctive features of 
solutions of density-matrix eq.uations for low and high 
strengths of irradiation, which have been fomd to be 
useful in the analysis, are discussed in Sec. IIG, along 
with a graphical method for interpretation of the spectra 
with low strengths of irradiation. 

Various mechanisms responsible for proton 
relaxation are described in Chapter III. It is often 
convenient to use the model of an isotropic external 
random field for the description of relaxation effects. 
This model includes the overall effect of various inter- 
actions, explicit considerations of which are rather 
involved. The conditions imder which various relaxation 
mechanisms can be approximated by a random field and the 
degree of correlation between these fields at the sites 
of different spins are discussed. It is shown that 
under conditions normally used in the experiments, 
several of the mechanisms can be approximated by a random 
field. Scalar coupling interaction is also examined and 
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tlie conditions mder wMcli it cazi be treated only as a 
mechazLism for relaxation are discussed. 

The freq_uency~ sweep homonuclear double resonance 
experiments were performed by adding auxiliary apparatus, 
to the Yarian HR-100 spectrometer with appropriate modi- 
fications. The details of the experimental set-up are 
given in Chapter lY. 

In Chapter Y, results of double resonance 

experiments performed on the ring proton (AB 2 ) spectrum 

of 2, 6-dibromoaniline (20^ solution in CCl^) are given. 

The parameters for the single resonance spectrum are 

^ ^ 89.95 ± 0.1 Hz. Double 

resonance spectra obtained for various strengths and 

frequencies of irradiation show significant relaxation 

featirces as well as marked variations in the inhomogeneity 

contribution to the various double resonance transitions. 

Two mechanism of relaxation (i) isotropic random fields 

and (ii) internal dipole-dipole interactions are considered 

for detailed analysis. The results of the analysis showed 

that the relaxation is primarily due to random fields 

having a ratio of mean square fields at ' A'' (position 4) 

and'* S'* (position 3 and 5) protons as 2:1. The random 

fields at the two ‘b** protons are found to be , 

almost completely correlated while those at proton ‘'A** 

c. ^ 

and any of the B protons practically uncorrelated. 
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Oiiapter VI oontains the proton relaxation studies 
in the AB system of 2-bromothi azole (neat liquid). The 
single resonance parameters ore = 3.55 + 0,05 Hz and 
I^A ” ^1 “ + 0,1 Hz, The scalar coupling of the 


protons with which has a fast quadrupole relaxation 
rate, gives rise to broadenirig in the spectrum and is an 
important source of relaxation of the protons in this 
molecule. The analysis of the double resonance spectra 
obtained for various frequencies and strengths of irra- 
diation leads to some significant information on the 


scalar coupling parameters including the magnitudes and 
relative signs of various proton nitrogen coupling 
constants, as well as the parameters describing the other 
mechanisms of relaxation in the sample. Conventional 
relaxation time and linewidth measurements were also 


used to facilitate the double resonance analysis. The 

relaxation is found to be primarily through scalar coupling 

and internal dipole-dipole interaction. The various 

parameters obtained by the double resonance analysis are 

for scalar coupling, sa- 9.7 Hz, Hz with 

AH sanie sign and for dipolar 

interaction t » 7.0 x lO"^^ sec. 
o 


The thesis has two appendixes, (i) a brief 
description and listing of a typical Fortran lY computer 
program and subroutines used in the density matrix calcula- 
tions and (ii) a reprint of an earlier work by the author on 
ME spectra of an ABXg spin system and solvent effects on 
proton- fluorine coupling constants. 
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I - liTG^EODUGOlIOlI 


lA SPIE EEI.AXATIOF-GENEML 

The interaction of a nucleus of magnetic moment 
ylfc. with a constant magnetic field H k, given hy^: 



' z o 


(I-l) 


causes a splitting of the (21+1) energy levels of the 
nucleus, corresponding to the eigenvalues of In a 

macroscopic sample^ consisting of an assembly of such 
nuclei, thermal equilibrium is reached by interactions of 
the nuclei with each other and with other degrees of 

»* M 

freedom in the sample, Icnovm as lattice , giving rise to 
a Boltzmann distribution of population amongst these 

levels and a net bulk magnetization along the direction 
of Hq. Transitions betv/een these levels can be induced 
by a radio-frequency (r.f.) field of appropriate frequency, 
applied perpendicular to The transition probabilities 

for upward and downward transitions are equal, and since 
there is an excess of population in the lower levels, this 
results in a net absorption of energy by the spin system 
from the r.f. field. This phenomenon is known as nuclear 
magnetic resonance. The absorption of energy tends to 
equalise the population of these levels and the rate of 


f The hamiltonians in this thesis are all expressed in 
angular frequency units. 



alDSorption is proportional to the population difference. 

On the other hand, the interaction of the spin system with 
the other degrees of freedom tends to restore the equilibrium 
population difference, She latter process is knovm as 
the relaxation process. Relaxation is manifest in the 
process by which a macroscopic sample developcsa bulk mag- 
netization = Xq the direction of the external 

A 

field H^k (where is the static susceptibility) and 
that by which the bulk magnetization in the x-y plane 
decays to zero, 

. A description of the relaxation process can be 
made by studying the time dependence of the components of 
bulk magnetization Ll(t), the equation of motion for which 
in the presence of a magnetic field H is given by Bloch ■’'s 
phenomenological equation [l]; 


dM 

dt 


= Y H - i 


^ _ i i 

^2 ^2 


h n 


( 1 - 2 ) 


According to this equation, for H = if H 

is in any arbitrary direction, the z-component of H appro- 
aches its equilibritam value exponentially with a time 

constant known as spin-lattice (or longitudinal) 
relaxation time and the x and y components of II approach 
their equilibrium value zero exponentially with a time 
constant '^ 2 * hnovm as spin-spin (or transverse) relaxation 
time. A determination of and Ig ^ts functions of external 
parameters can be used to derive information on the 
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interactions responsilole for relaxation. In complicated 
systems, like liquid samples giving cmltipX-e line epectra 
Tindex Mgli methods like double resonance are 
used for study of relaxation processes. The present v/ork 
concerns itself with a study of relaxation processes in 
coupled spin systems in liquid phase by steady-state double 
resonance. Before describing the double resonance techn- 
ique, the particular advantage of the method for the 
study of relaxation in coupled spin systems in liquids is 
briefly outlined in the following. 

The characteristic features and the methods of 
experimental study of relaxation are different for different 
situations. In solids the resonance lines are usually 
broad due to dipole-dipole interactions, leading to a 
short spin-spin relaxation time T2. The spin-lattice 
relaxation time T^ is comparatively long. The spin-spin 
relaxation may be studied experimentally through linewidths 
and lineshapes of resonances v/hile, in general, both T^ 
and T2 can be measured by experiments in which the transient 
response of the spin system to an r.f. perturbation is 
monitored (e.g. spin-echo method) to determine the relaxa- 
tion times on the basis of Eq. (l- 2 ). 

The situation is rather different, by contrast, 
in liquids. In liquids the dipole-dipole interactions are 
averaged out by the rapid molecular motions leading to 
extremely sharp resonance lines, the widths of which are 
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primarily determined by tlie inliomogeneity in tiie applied 
magnetic field Furthermore^ the sharpness of the resonance 
lines in liquids enable the detection of small variations 
in the electronic environment (chemical shift) and interac- 
tions arising from indirect spin-spin coupling in a molecule, 
which lead to a fine structure in the high resolution n*m.r. 
spectra. The frequencies and the intensities in such 
spectra can be described by the hamiltoniaHj for a single 
molecule^ given by [2]: 

( 1 - 3 ) 

where v^^ “ “ "Yi \7hich the gyromagnetic 

ratio including the chemical shift of the i*« nucleus,' 
and J . . is its spin-spin coupling constant v;ith nucleus 3. 

X J 

Transient experiments like the spin-echo method can be used 
for measurement of relaxation times in liquids if the 
spectra flo not exhibit multiple t structure. TJhen there 
is multiple t structure, these measurements encounter 


f Exceptions to this statement occur for the resonances 
of nuclei with spin ^ 1, relaxing due to quadrupolar 
interaction^, or for resonances of spin ^ nuclei coupled 
to either quadrupolar nuclei or nuclei participating 
in chemical exchange etc. 
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several difficulties. Separate T2_ -and I2 measurements on 
eacii line in a multiplet are not usually possible because 
the lines are closely spaced, and only an overall effect 
of a number of lines may be obtained. Furtheimore, the 
description of the relaxation process by single exponen- 
tials with time constants and !]}2> Bloch equations, 

is inadequate. Therefore, for a study of relaxation 
processes in such cases, one needs a general theory and 
an experimental method, the results of which are dependent 
on relaxation processes. T/angsness and Bloch [ 3 ] and 
Bloch [4] have given a density-matrix theoretical descri- 
ption of the phenomenon of nuclear induction. The density- 
matrix treatment is quite general and suitable for 
description of relaxation effects in a variety of systems. 
Double resonance experiments, the results of which are 
interpreted by the density -matrix theory, meet the above 
requirements for relaxation studies in coupled spin systems 
in liquids. 

IB DOUBLE EESOUANCE 

The double-resonance experiments consist of 
irradiation of one of the transitions in the spectrum by 
a strong r.f, field and simultaneous observation of the 
remaining parts of the spectrum by another weak r.f. field. 
The principle of this method is based on the fact that the 
strong r.f. field competes with the processes of relaxation 
in the spin system, thereby causing a variety of effects 

viz. changes in populations and lifetimes of the energy 
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levels which in turn are reflected as changes in intensities 
and linewidths in the double -resonance spectrum [4]. The 
quantitative details of these features depend on the rela- 
xation processes in the spin system. Information on the 
relaxation processes is then obtained by comparing the 
observed features with those calculated for various rela- 
xation mechanisms. 

Strong irradiation also mixes various eigenstates 
of the spin hamiltonian (Eq. 1-5) , giving rise to new 
transitions and shifts in frequencies of single resonance 
transitions [4]. These features are implicitly included 
in a detailed analysis of double -resonance spectra even 
though they do not depend on relaxation, and can be 
completely undei stood by a diagonalization of the ^in 
hamiltonian including the irradiating field, in a frame 
rotating at the irradiating frequency [5^1. 

A density -matrix description of the phenomenon 
of nuclear induction as formulated by Bloch [4] and 
Redfield [6], and given by Abragam [7] and Hubbard [8], 
is capable of describing all the effects of the double 
resonance experiments including, in particular, the 
relaxation effects. In this formalism the dynamical 
behaviour of the spin system is described by rewriting 
the general equation of motion of the density matrix, 
using second order perturbation theory for the spin-lattice 
interaction. This results in differential equations for 
the various elements of the density matrix, which include 
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the effects of relaxation through the elements of a general 
'‘relaxation matrix^'" [6] . The various relaxation matrix 
elements contain, for a multilevel system, the information 
implied in and Tg Bloch’s phenomenological equation 
(1-2). The results of double-resonance experiments are 
then described through changes produced in the density 
matrix due to irradiation. This involves a consideration 
of the complete relaxation matrix. A calculation of these 
changes for various probable relaxation mechanisms is then 

i 

compared with the observed features in a steady-state 
double -resonance experiment. This method thus avoids 
some of the theoretical ambiguities and experimental 
difficulties in and T2 measurements, mentioned earlier. 
Furthermore^ these experiments provide an example in which 
a detailed consideration of the density matrix formalism 
is required and thus form an experimental test for the 
theory. 

In the following some of the previous vfork on 
the study of relaxation in coupled spin systems by the 
double-resonance technique is briefly described in order 
to discuss the motivation and the nature of the work 
presented in this thesis, 

Baldeschvmeler [5] gave a detailed method of 
calculation of relaxation effects in double resonance 
based on the density-matrix formalism of Bloch [4] and 
Redfield [6], This method was applied to the proton double 
resonance in the single spin system of CHCl^ and some 
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discrepancies were noticed between experiment and theory. 
This method of analysis was later on re-examined by Hageswara 
Eao [9] and a similar but somewhat improved method was 
suggested. It was shovm that the approximation originally 
proposed by Bloch [4], for systems with weak relaxation, 
for the analysis of the double resonance spectra obtained 
under conditions of strong irradiation, is applicable 
even when the frequency of irradiation is very close to 
a single resonance transition. A reinvestigation of proton 
double resonance in CHCl^, on this basis, removed the 
discrepancies between the theory and the experiment for 
the single spin -1; case [10]. 

Nageswara Rao, Baldeschwieler and Anderson [11], 
used the above method to analyse the results of steady- 
state frequency-sweep double resonance experiments on a 
strongly coupled two spin system (AB) fomed by the protons 
in 2-bromo-5-chloro thiophene, under conditions of strong 
irradiation. The validity of Bloch approximation for 
strong irradiation spectra has been verified in this 
analysis and furthermore it was possible to clearly dis- 
tinguish between different relaxation mechanisms on the 
basis of a comparison between the observed and' the calcula- 
ted relative intensities in the double resonance spectra. 

The relaxation in this molecule was found to be predomi- 
nantly due to correlated external random fields. Wakefield 
and Memory [12] recently made similar studies on an AB and 
an ABX type of spin system, using Bloch approximation, 
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and concluded tliat the correlated random fields fit their 
spectra best. Steady-state double resonance experiments 
have also been performed on an AB spin system by Matsuoka 
and Kakiuchi [ij] , to obtain quantitative information on 
the contribution to relaxation for various intramolecular 
and intermolecular interactions and concluded that a single 
exponential cannot be used to describe various relaxation 
rates. 

Hoggle performed transient as vfell as steady-state 
double resonance experiments on the AB spin system of 2,5- 
dibromothiophene [14] . Ihe transient experiments involved 
a periodic observation of recovery of single resonance 
intensity after the saturating r.f, field is turned off. 

He also examined the correlations between external random 
fields and concluded that both internal dipole-dipole 
interaction and external random fields with partial corre- 
lation are responsible for relaxation in this molecule. 
ICuhlmann and Baldeschwieler [15] made similar studies in 
difluoroethylenes. Ihese authors used low strengths of 
irradiation such that line splittings and shifts do not 
occur and only intensity changes in the single resonance 
spectra are observed (generalised Overhauser effect [4]) * 

Ihe spectra v/ere analysed by using external random field 

model for various relaxation mechanisms. It was concluded 
that spin-rotation interaction is responsible for relaxa- 
tion between states of different •.symmetry. In this case 
information regarding the relaxation rate of fluorine 
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could Tse obtained v/liile observing only proton resonance. 

Steady-state double resonance studies have been 
extended by Hageswara Eao and Lessinger [16], to weakly 
coupled systems containing a large number of spins. It 
was shorn that if the interactions causing relaxation in 
such systems satisfy certain symmetry requirements, it is 
possible to obtain quantitative information on the relaxa- 
tion by assTUiming a random field model. The double reso- 
nance experiments were performed in etliylfluoride and the 
results indicated that the relaxation of fluorine is 
approximately four times as efficient as that of protons 
in this molecule. 

Closely connected with the studies mentioned above 
are some hetronuclear double resonance experiments in the 
AX spin systems of G in chloroform [17] and that of 

H -D in gaseous HD [18] . Significant relaxation infor- 
mation has been obtained by using the density-matrix 
analysis for these systems. Interesting extensions of 
the steady-state and transient double resonance experi- 
ments have been made by Gordon and Baldeschwieler [19,20] 
and Gordon [21] in which the observing and the irradiating 
r.f, fields are periodically pulsed. Double resonance 
experiments have also been performed in liquid mixtures 
(inte molecular Overhauser effect) [22]. 

These experiments clearly demonstrate that the 
double resonance experiments coupled r/ith density-matrix 
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analysis are capable of yielding wortliwliile information 
on the relaxation processes in different types of spin 
systems and under a variety of experimental conditions. 

IC PEESENT WOEIC 

The work presented in this thesis, on coupled 
spin systems, is undertaken to study the systematics of 
the double resonance method and to obtain quantitative 
information on the various mechanisms of proton relaxation. 
The scope of the present work is restricted to steady- 
state double resonance experiments. The work is based on 
the study of proton-proton double resonance spectra of 
(1) AB 2 spin system formed by the ring protons of 2,6- 
dibromoaniline studied as a 20 per cent solution in CCl^ 
and (2) AB spin system formed by the protons in 2-bromo- 
thiazole studied as a neat liquid. Double resonance 
spectra wexe obtained, in each case, by irradiating each 
one of the single resonance transitions over a wide range 
of strengths of irradiation. 

Prom the analysis of the AJB 2 spin system the 
following points have emerged: it is found that an explicit 
consideration of the effect of ■ inhomogeneity of the 
magnetic field is of considerable importance in the 
determination of quantitative information on relaxation 
from double resonance spectra. Freeman and Anderson [23] 
have shown that the inhomogeneity contributes, in general, 
different widths to different transitions in a double 
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resonance spectrnn. In -the present worls: a method has "been 
suggested to include this in the density matrix analysis 
of double resonance spectra and this led to a reliable 

* 

estimate of the relaxation parameters. This effect has 
been overlooked in earlier studies of relaxation by double 
resonance [11-14 ,1S yl8, 24] . Secondly^ it is shown in this 
work that it is often useful to study the double resonance 
spectra over a wide range of strengths of irradiation. 

This is due to the fact that the solutions of the density- 
matrix equations for high strengths of irradiation are 
determined by the relative values of relaxation matrix 
elements while those for lov/ strengths of irradiation 
depend on their absolute values. This distinction is of 
considerable value in distinguishing between different 
mechanisms. 

IEt 70 relaxation processes (i) external isotropic 
random fields and (ii) internal dipole-dipole interactions, 
along with various probable degrees of correlation in each 
case, were considered for a detailed analysis of relaxation 
effects in the ABg spin system, The correlation constants 
for the isotropic random fields are defined to include the 
possibility of complete correlation accompanied by different 
mean square fields. The analysis led to the conclusion that 
the relaxation in the spin system is primarily due to 
external random fields with a ratio 2:1 for the mean square 
fields at the ''' k''' proton (position 4) to those at the^B'"'' 
protons (positions 3 and 5) , The random fields at the two 



14 


B proton sites are fonnd to Be highly correlated while 
those at the ^a' proton site and any of the ‘b*^ proton sites 
are completely nncorr elated,"^ 

In the AB spin system of 2-hromothi azole scalar 
coupling of protons vd-th which has a fast quadrupole 

relaxation, was found to modify the process of relaxation of 
the former to a large extent. This coupling causes significant 
Broadening in the proton spectrum of this molecule. The 
douBle resonance study yielded information on the various 
parameters descriBing the scala.r coupling interactions viz,, 
the magnitudes and relative sign of the proton-U^^ coupling 
constants. The scalar coupling interaction was also 
examined in some detail and the conditions under which it can 
Be treated only as a mechanism for relaxation are discussed. 

It is found that under the conditions of the experiment the 
scalar coupling interaction contriButes significantly only 
to spin- spin relaxation, so that the contact of the spin 
system with the lattice should Be through other interactions. 

A detailed douBle resonance analysis for Both high and low 
amplitudes of irradiation cases lod to the conclusion that 
the internal dipole- dipole interactions along with tho 
scalar coupling account for the spin relaxation in this 
molecule. A graphical method for interpretation of the 


-j- The results of this study are puBlished in a paper By 
Anil Kumar and B.B. Bageswara Rao, Molec. Phys., 15 , 
377 (1968). 
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V t 

spectra with low amplitudes of irradiation is adopted for 
this analysis.'^ 

It is often convenient to use the model of an 
isotropic random field for the analysis of relaxation 
effects in the double resonance studies. This is because 
this interaction is linear in the spin operators I_, I_ and 

X y 

I and the lattice variables and possesses isotropy. This 
z 

model represents the overall affect of various intramolecular 
and intermole cular interactions explicit consideration of 
v/hioh is rather involved. As part of this work, the condi- 
tions under which different mechanisms can be approximated 
by a random- field model are examined. It is shown that 
under conditions normally used in experiments, several of 
the mechanisms can be approximated by such a model. 

The experiments wero performed at room temperature 
and by frequency- sweep method. The experimental arrangement 
is similar to the standarad arrangements used in literature 
[23,25] with some minor modifications arising from the 
particular apparatus used. This arrangement was set-up and 
was used along with a Varian HR -100 n.m.r. spectrometer. 

The plan of presentation of the thesis is as 
follows. The notation and the density-matrix method used 
in this work are given in Sec. IIA. The treatment of 
inhomogeneity and its contribution to single and double 

^ The results of this study form part of a paper by Anil 

Kumar, H. Rama Krishna and B*D. Hagoswara Rao, Molec. Phys. , 
1969 (in press) , 
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resonance spectra are given in Sec. IIB. The solutions of 
the density-matrix eq_uations for different strengths of 
irradiation are discussed in Sec, IIG along with a graphical 
method for interpretation of spectra for 'low^ strengths of 
irradiation. Various relaxation mechanisms are listed in 
Sec. IIIA, and the conditions under which they can he 
approximated by a r-andom field are examined in Sec. IIIB. 
Chapter III also contains a discussion of the conditions 
under which scalar coupling can be treated only as a mechanism 
of relaxation. The details of tho experimental arrangement 
for frequency- sweep double- resonance experiments are given 
in Chapter IV. Chapter V contains the results of double 
resonance experiments on bhe AB 2 spin system of 2,6-dibromo- 
aniline and their analysis. Double resonance spectra in 
this molecule showed significant relaxation features and 
marked variations in the contribution of inhomogeneity to the 
linewidths. Those features are discussed in detail. Chapter 
VI contains tho results of the double-resonance experiments 
and analysis of the relaxation effects in the AB spin system 
of 2-bromothi azole. The conclusions and results of the work 
are summarized in tho end. There are two appendixes to the 
thesis (i) a brief discussion and a listing of a typical 
computer program used in the density-matrix analysis and (ii) 
a reprint copy of an earlier work on solvent effects on 
coupling constants, by the author. 
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II - THEORT 


IIA BASIC DEBSIIY MATRIX EOEMAIISM 


1. Equation of Motion 

The density-matrix method of analysis of double 
resonance spectra given recently by Baldeschwieler [5] 
and Rageswara Rao [9l> based on the formalism of Bloch [4] 
and Redfield [6], (also discussed by Abragam [7])j is 
discussed briefly in this section. This presents the 
basic assumptions involved in the theory and the notation 
to be used in the later sections. 

The equation of motion of the density matrix 
of the total system (including the nuclear spin system, 
the lattice and thoir interaction)^ in Schrbdinger repre- 
sentation^ is given by; 


df 

dt 


0 nf 




(Il-l) 


where is the spin hamiltonian which is time dependent 

because of the observing and irradiating r.f. fields, 
is the hamiltonian of the lattice and is the spin- 
lattice interaction. 

Transforming Bq. (Il-l) into an interaction 
representation by defining ; 


= U(t) e^^^ f U“^(t) 

^'^(t) = U(t) e^^'^ e-iP't u“^(t) 


(II-2) 
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where U(t) = U(t,0) obeys .the differential equation 


= i IJ(t) ^(t) (II-3) 

gives: 


a?^ 

at 




(II-4) 


Equation (II-4) solved in two successive approximations 
leads tos 


0 

J^(0)]]dT, (II-5) 

A spin density matrix cr is defined hys 

a^(t) = Ir^ ( ^'^(t)] , (II-6) 

where 

a^(t) = U(t) a(t) U"^(t) (II-7) 

and \ represents the trace over lattice states. 

Equation (II-5) is solved with the following assumptions: 
i) The time t = 0 is so chosen that the system is in 
thermal equilibrium at that time. The lattice is assumed 
to have a large heat capacity and to remain always in 
thermal equilibrium, ^ (O) is then factorizable into spin 
and lattice parts: 
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f^(0) = >J(0) cr^(O), (II-8) 


where yj(0) = \(0) and cr^(O) = cy(0) are the equilibritim 
density-matrices at the lattice temperature for the lattice 


and spin parts respectively. 


^ e 


, known as relaxation hamil> 


tonian, is a random operator uncorrelated with a(0). Ihe 


trace over the lattice states will then give zero for the 
first term on E.II.B. of Eq.. (II-5). 

iii) The upper limit of integration in Eg, (II-5) can he 


extended to infinity if t » v , where t is the correlation 

'' c c 

time for time dependence of ^(t), T/liile t » t it 
should, however, he short enough to satisfy the perturbation 
expansion of Eg, (II-5). Eor short correlation^ a^(0) in 
the second term on the R,H.S, of Eg.(II“5) may he replaced 
by cr^(t) , 

iv) The time independent part of is such that 

it ^^/kT 1 (high temperature- approximation). This allows 

T T 

a replacement of o (t) by (o (t) - o(0)) in the second 
term on the R.H.B, of Eg. (II-5), while performing the 
operation ^ 


Transforming back to Schrodinger representation 
and noting that in liguids^ where the l^^ttice motion may 
by adeguately described classically, the operation, of 
trace after multiplying v/ith JX(0) can be replaced by a 
classical ensemble average [26], denoted by ^^^one 
obtains?, 
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= -if 
dt 


where 


(_(t) j a] - 


ft*’ 


t) , [U(i:“T:,t) 



t-T) 


ir~^(t-T 5 t), 0 ( 1 ;) - cr(0)]] ' (11-9) 


cf(0) - 0^ _ .J; (1 ^ 




-- 1 

1 


•) = 


q. 




(II-IO) 


1 is the niimher of spin states and q. = !fc./lTMl, 

At this point one of the following two assumption 

is usually made, depending on how short t is; 

o 

(a) 'Extreme narrowing This assumption recquires 




o^c 


«1. If it is valid_,U(t-T,t) of Eq.. (II-9) 
becomes unity and the spectral densities (to he defined 
later) become independent of freq.uency [ref. 7, p. 280]. 

(b) Eon-viscous liq.uid'*' approximation. If t is not 

c 

short enough to satisfy (a), it is most lihely that 
I ..4 ^^ 2 ^(t) <K 1, where the summation k is over all the 


k 


time- dependent parts of spin hamiltonian (In double 

1Kcse«.ve cLuue.i:o 

resonance^ the observing and irradiating r.f. fields, 
end In all the cases of present interest^ the strengths 
of these fields are very much smaller than [ref. 7, 

p. 515 ]). Ihis leads to U(t-T,t) = exp (- in 

Eq. (II-9). 

Although the exact form for depends upon 

the particular interaction between the spin system and 
the lattice, it can^in general^ be expressed as: 
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‘ir 

wliere A^'^^ are spin operators and I'^‘'-^^(t) are random 
functions of lattice variables expressed in irreducible 
form. The liermiticity of requires s 



(t) 


(_ 1)1 and A^^^^ = (-1)^^ (II~12) 


This form for 



\ 



to the simplification: 



\ 



(11-13) 


known as '•'secular approximation'^ , vdiich arises from the 
fact that the ''non-secular" terms (q q-' ) vary rapidly 
with time and give zero on averaging over the ensemble[ 4 ]. 
Using Eqs. (II-ll) - (II-I 3 ) and the *■' non-vise ous liquid’"'^ 
approximation, Sq. (II-9) becomes [25]: 


= -i[ ^^(t), a] - :^(t), [ ^^(t-x), a-c^]] 


^^ 0 - • _ jJi: 

dt 


dlt. 

/av ^ 


(11-14) 


where 


3t(t-T) = 


(-I)*! , 


(11-13) 


in which 


Y 

i 2 tt;xv^ V I^( j) 

® > 

(11-16) 



25 


■where and in writing Eq, (11-16) the coupling 

constants and the chemical shifts have been neglected 
compared to Xarmor frequencies in '■'eztreme- 

narrqwing'^ limit the equation of motion of cr, Eq. (11-9), 
becomes [ 9] 5 


_ 

dt - 



CT] 



(11-17) 


2. Double Resonance - Rotating Coordinate Iransf ormation 


Ihe spin hamiltonian; 


consists of a time- independent part defined in Eq.(l-3), 

and interactions of the spins with a wealc observing r.f. 
field (/l^^(t) and a strong irradiating r.f. field f 

given by; 


^l?j^(t) = D^.^ esp(icO^^t) + Dj^_ exp(-iC 0 i_t) ^ (11-19) 


where 


and 


V. 


ki 


= - YA/2^t 


k =1,2 


( 11 - 20 ) 


( 11 - 21 ) 


in which and v^.^ are expressed in Hz, and ^ = H^cos 
-H^ singj^t 3 is the observing r.f, field and H 2 = 

A A 

cos (j gt i - Eg sin 6 ^ 2 '^ 3 irradiating r.f. field. 
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T rails foxiiiinr to a coordinate frame rotating at the 

/V 

freq.uency (-(hg];) of the irradiating field H 2 ihrongh 
the operator 


1 


exp ( -iol 2 't 



and noting that: 


( 11 - 22 ) 



^^(t) = exp(ioo't) + I)^_ exp(-ico^t) , 

[ f [“o’ 5^] = 0 ; 


( 11 - 23 ) 


where denotes operators in the rotating frame and 
UJ = hq. ( 11-14) gives: 


dt 


= -i[^? + - _r<[ ifft), [ 


\ 


where 


JL 


R 


/ XJ 

1 — 1 

1 — 1 

0 

b 

(11-24) 

^^(t), n] ~ [7? - o^) 

, (11-25) 

^ 12(1) + ^2 


+ Z7' Jii I(i)«I(3)] t 

i < D ^ 



(11-26) 

= + (^2/2.1^, 

(11-27,) 


in which 
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(a - a^) is. the relaxation term defined Tdj its matrix 
elements in a basis in wMch is diagonal if tlie solution 

is obtained in tlie rotating frame tiirougii Eq.. (11-24) . 


(cr - a^) will similarly be the relaxation term defined by 
its matrix elements in a basis in which is diagonal if 

the solution is obtained in the laboratory frame through 
Eq. (11-14). In general in a basis a, a'.... 


V ' 

\ 


aa 


where 


,(a - a^) = > 


(11-28) 


r"-*— > 

^aa'' ~ “ ^a.B ^ — ‘ '^Ya''v6' "" ^ 


in which 


^apa''p^ 


apa'p 


p'a'pa ^ap P'4^ " "* > 

(11-29) ' 


jAj ^(t)!p) <^a'! ^ 


At. 


av 
(11-30) 


In writing the spectral densities this form, 
the imaginary part, which is usually small, is neglected 
[ref. 7, p. 279]t. 

The hermiticity of #t) yields [9]j'^’*’ 


f This term is zero for the “extreme-narrowing*^ assimption. 

H 

These symmetry relations are obtained for extreme narrowing 

XV 

and will hold for non-viscous liquid approximation 
if the chemico.1 s'aifts and coupling constants are negle- 
cted compared to jLarmor frequencies in writing the Eq. (11-16). 
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~ ^p^a^pa » 

giving; 

\a''pp^ " \^ap'p ~ ^pp'aa^ “ ^p'pa^a ' 

and 

^ \app' ^ 

Cu 


(11-31) 


( 11 - 32 ) 


(11-33) 


Using Eq.s. II“(ll-13), (11-15) and (II-16) \?e get for 
Eq.. (II~30) under'^non-viscous liq.uid^'" approximation; 


= 1/ (ajA(-4)jp> 

q -oo . / 


(v) I P^ 


(11-34) 


and for '''"extreme -narrowing'^" approximation: 


rsn _ p 

'^apa'p' " 


<1 


E 


(l), 


“> av <“| 


I (-a)i \ 


A' 






(-< 1 ) 

(11-35) 


3 . Solution of Equation of Motion 

Equation (11-25) is solved in tvro steps, first 
by determining tUe deviation ”X(’'^) caused by the irradiating 
r,f, field ^^^ 2 ^”*^) unperturbed siain density matrix 

a and then by obtaining the deviation u(t) caused by 
the observing r.f, field on t]ie resulting system. 

The total spin density matrix, including the effect of two 
r.f, fielriSjis then expressed as; 
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cr(t) = Oq + X(t) + riCt) . (11-36) 

In writing this expression, linear superposition of the 

« 

effects of the irradiating and the observing fields on the 
density matrix is assumed. 

In the rotating coordinate systemi 


'^(t) = cJq + X + ^C'fc)/ (11-37) 

where X is independent of time, since ^^2 independent 
of time in the rotating frame. 

Effect of Irradi ation s- Setting ^^!^(t) =0 ; 'n(t)=0, 
in the first step of j^solution, Eq. (11-25) becomes in the 
steady-state ; 


[ 





( 11 - 38 ) 


Using Eqs. (II-IO) and (11-26) and neglecting the chemical 
shift between irradiated nuclei compared to their --Carmor 
frequency one gets on taking the matrix elements of Eq.(lI-38) 
in the a,a^.... basis ( such that | =E^|c^ ): 


(E 


a 




+ i 


i: 

pp" 


R 


aa' 




■{2%q,v^) (uYg) 


[1^(3) 3 ) laa'' * (11-39) 

D 


Since the densitj?- matrix is hermitian Ee('Y y) = Re('X r ) 

^aa. ' ^ a ct^ 

and Im (X^.^^/) = -Im an^ ihe diagonal elements of 
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% are real* Separating the real and imaginary parts of 
Eq. (11-39) gives for a ^ 




l^'< P 




-PP' 


= -(STxqv^) (TCVg) Y1 [I+(d) - (11-40) 

«} 


^^n'.n'/RR’' ■*' ^’n-n' '’t/r^ R®( 'XfiR''') 


p<p 




^PP 


• I I ■ «. 


(11-41) 


and^ for a = a'^; 

^ \app' ^ \app Xpp = (11-42) 


Equations ( 11-40) to (11-42) form a set of simultaneous 
equations in the matrix elements of v;hich the I 

equations given by Eq, (11-42) are linearly dependent. 
Using the condition Ir a = Ir cr^ = ly gives: 



(11-43) 


and a unique solu oion for the complete 3^ matrix can be 
obtained by replacing one of the equation in ( 11-42) by 
Eq. (11-43). Eurtlier simplification of the set of Eqs. 
(11-40) to (11-43) is possible depending on the strength 
of irradiation. I’his aspect v/ill be discussed in Sec. IIC. 
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Effect o f OlDs ervi ng Field :- Ihe deviation in tlxe spin 
density matrix (a^ + ^ ) • due to observing r.f« field '^^(t) , 

/'V*' 

given Toy il(t), is obtained using Eqs. (11-25) > (II-37) and 
(11-38) as: 


dp(t) 

dt 



^(t)] -i [#,(t) 


+ 7- + E(t)] 


(TT(t) ) , 


(11-44) 


A steady-state solution for T](t) for small values of H^(t), 
small enougli not to cause any saturation effects, is obtained 
by assuming that the diagonal elements of 'n(t) are zero and 
that the off-diagonal elements have a time dependence given by: 


Tl / 

aa 


(t) 


„+ ^i t - 
aa aa 


-i t 


(11-45) 


Using the hermiticity for p(t), substituting (11-45) in 
Eq.. (11-44), collecting terms varying as esp (i<^^t), 
exp (-i'-^'^t) and constant in time, and separating the real 
and imaginary parts give a set of linear simultaneous equa- 
tions in the matrix elements of 'n(t). If is not very 
close to small H^(t), these equations simplify 

to give [9]: 


Ee(Y±„,) + ) Im(Y± = 


p<r 


-ChH- ’ 


( 11 - 46 ) 



i 


30 


-f 


p<p-' 


■[“l+ ’ ^ ]aa' , 


( 11 - 47 ) 


i 


in wiiich tlie upper sign is important if iu>^+ is 

small and the lovrer sign if (- caj''+ E -E / ) is small. 

CX vJCi 

Only one of the off-diagonal elements of Ti(t), 
corresponding to the transition "being observed, is important 
at a time if there are no overlapping transitions present 
in the double resonance spectrum, 0>\€ likeivob tains: 






aa'" 


Ea - + Ka'aa' 


(II-+8) 




-\a^aa' Chi> (E® X+%) -(±cjVe„-B„.) [Bl^, Im 

(±W'+ Ea - V)' + C-acc' 


( 11 - 49 ) 


^ • Signal Intensity 

The macroscopic expectation value of an operator 
is given by a trace of the product of the operator with 
the density matrix, and therefore the signal in a nuclear 
induction experiment is given by [5,9]: 
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A T 

'21 dt 


Yi 


i i^pt 

e 


O' 


I (i) 


i u)_t 
e 


a 


} 


( 11 - 50 ) 


where K^is a proportionality constant. The diagonal 
elements of a which are time independent contribute signals 
at aip, and the off -.diagonal elements at (+ » <1^, 

I ^ c X 

O-vid 

(-W^+ 2 i-v 2 '* intensity of a transition at i>'>^,in 

the absorption m.ode^ is given by: 


i 

= +■ ^ Tj_ sin co^t 

i 


"^'acc ^^1+’ ^'*'‘^ 0 ^ ^aa' ^ 

^ ^ r'- . ..... .r':: j 

(±6.''+ \ + RL/../ 


aa 'aa" 


(11-52) 


The second terra in the numerator of Eq. (II-52) is negligible 
if Bloch appronination holds (to be discussed in Sec IIG), 

■i- 

and the signal S~'^, / becomes a lorentzian with half-vridth 

i -n 

\a'^aa'' hali-height. The peak signal intensity is 

given by: 
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S- 

aa 


± K( 4-Yi(I+{ i))a<x/ + X+%) ]^a' 


■p 

aa’^aa'' 


(11-53) 


where a redefined K(tO^,H^) includes all proportionality 
constants. The integrated intensity (integrated over 
AUJ) from Eq. (il-52) is obtained as; 


Taa' = ± K('Ji.Hi) 






= (llnewidth)^^, - 


(11-54) 


EquAuon (11-54) is correct eTen if the second terra in 
the numerator of Eq. (II-52) involving Im(X) is not 
neglected, since the corresponding integral vanishes. 

This term thus superposes an intensity distribution odd 
with respect to A 60 on the lorentzian, but leaves the 
integrated intensity unchanged. 

The expression for peak-height jEq. (11-53), 
contains the relaxation width in denomi- 

nator, but in a high resolution experiment the lineshape 
is determined by a few other factors, such as the 
inhomogeneity of the magnetic field which contributes 
significantly to linewidths. These factors should be 
considered before making a comparison of the peak-heights 
with the observed spectra. This point v;as not explicitly 
considered in earlier studies of relaxation effects by 
double resonance [11-14,16,18,24]. 
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IIB MAGNETIC EIELD INHOMOGENEITY 
1. Su’b-ensem'ble Treat ment of Density-matrix 

If the magnetic field is inhomogeneous over the 
sample volume, the spin Hamiltonian is different for 
different spin-Isystems in the sample and the treatment of 
the previous section should then he modified to include 
this. Ahragam [ref. 7, p. 26] suggested a method of 
handling this problem by constructing sub-ensembles, such 
that the spin-hamiltonian is the same for the systems in 
each of them, but varies from one sub-ensemble to another. 
The number of systems in each sub-ensemble should be large 
enough to treat the problem statistically, but small 
enough, at the saae time to neglect variations in the 
spin-hamiltonian vmthin a sub-ensemble. The density- 
matrix theories given earlier [4-9] and sketched in Sec. 
IIA here, will then be applicable to each sub-ensemble and 
a final result can be obtained by appropriate averaging 
over all the sub-ensembles, 

For each sub-ensemble one can then define a^ 
density-matrix , the equation of motion of which is 
given by; 

, (11-55) 

where is total hamiltonian for the sub- 

ensemble and is^same for all the spin-systems within it. 

The time variation of is then given by: 

J 
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J’j = U.(i:) ^.(0) UjHt) (11-56) 

where U.(t) satisfies the differential equation (see Eq, 

J 

II-3) : 

dU.(t) • nn 

An important point to note here is that if the density- 
matrix y for the entire ensemble is defined by averaging 
over the density-matrices of individual sub-ensembles, 
the transformation connecting ^(t) and ^(0) may not, in 
general, be unitary, even though that connecting 
and 5d^0) each sub-ensemble is unitary [ref. 7, p.27]. 
However in the present problem it is not necessary to 
define an average density matrix and the entire double 
resonance calculation may be performed for each sub- 
ensemble by solving the equation of motion of y., given 
by Eq. (11-55) , using the theory sketched in Sec IIA, The 
relevant macroscopic quantities, like the time derivatives 
of magnetization^ can also be computed for each sub-ensemble, 
and a final result may be obtained by averaging these quan- 
tities over all the sub-ensembles. If the magnetic field 
inhomogeneity may be assumed to satisfy the requirements 
of such a treatment then the signal intensity for the 
total ensemble is obtained by integrating the intensity 
obtained for a sub-ensemble given by Eq. (II— 52) , after 
multiplying it by the shape function of the inhomogeneity. 
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Assuming that the shape fimction (normalized) 
for magnetic field inhomogeneity h(A is a lorentzian, 
given hy: 


h(A O) 


l/itT^ 

(AfA3)2 + (1/2*)^ ^ 


( 11 - 58 ) 


where l/T^ = YAHqj half-width at half-intensity 

due to field inhomogeneity, multiplying it with Eq_. (11-52) 
and integrating over A f 2 ?om -COto +oo , leads to a 
folding of two lorentzians, which gives another lorentzian 
having a half-v/idth equal to the sum of half-widths of 
each component [ref« 7, p. 50]. 

In addition to the relaxation width and the 
inhomogeneity width, there is an extra linewidth arising 
from the experimental set up, due to factors like minute 
instabilities in the audio-oscillators [23]. This width^ 
though small, is an unavoidable factor, v/hich may be 
assumed to be the same for all the transitions and may 
be added to the linewidth, in the expression for intensitj^ 
in an adhoc manner. The peak-height of a non- overlapping 
transition in a double resonance spectrum is then obtained 
as: 





(11-59) 
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where is the total half-width of the a ->> 

transition given hy; 

~ ^aa''aa'' ^ Inhomogeneity '^^^Insta'biliiy 

(11-60) 

The peah-height in the single resonance spectrum 
is then given by; 

S- i = , (11-61) 

where a, a''.... is the basis in which is diagonal and 

consists of Haa^aa-’ ^ At^aa- )l^omogenelty' 

( AW) ^ 

2. Inhomogeneity Contribution to Linewidths 

The linewidths in high resolution n.m.r. spectra 
of weakly relaxing spin systems are primarily determined 
by inhomogeneity in the magnetic field, Por single 
resonance spectra this leads to the situation that all 
the transitions have approximately the same observed 
linewidths inspite of any variation in relaxation widths 
arising from the nature of the particular relaxation 
mechanism for the spins in the molecule. Por double 
resonance spectra, however, the contribution of the 
magnetic field inhomogeneity oan be different for 
different transitions as shown by Preeman and Jhiderson [25].; 
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This is basically due to the fact that the process of irra- 
diation of one of the nuclei in a molecule^ with a monochro- 
matic freq[uenGy, causes freq,uency shifts of either sign in 
the transitions of other coupled nuclei, the magnitudes of 
these shifts depending on the frequency offset** of 

* ^ vjt 

irradiation. The freq.uency offset depends on the local 
field at the irradiated nucleus due to the ihhomogeneity 
of Hq, and these local fields at both the observed and 
the irradiated nuclei are exactly correlated [27], In an 
inhomogeneous field the '''' frequency offset*'’' has a distribu- 
tion of values giving rise to a spread in the shifts of 
the line positions, due to irradiation, which may either 
cancel or augment the effect of inhomogeneity, giving 
rise to lines which are narrovrer or broader compared to 
s ingl e res onanc e . 


(a) Single Eesonance :- In general, the width caused 

by an inhomogeneous field of width in a transition 

a a'*', may be written as: 


C 


Ihbomogeneity 



± * 


(11-62) 


' 1 ' 

where is the frequency of the a -> transition and 

Cu^ the Xarmor frequencies of various spins in a molecule. 
For a weakly coupled single resonance spectrum Eq. (11-62) 
gives III "the case of strongly coupled single 

resonance spectrum, the farmer frequencies of other nuclei 
enter the expression for through the chemical 



shift (6j- - Lo-) and if all the strongly coupled nuclei 

1 J 

belong to the same species, Eq.. (11-62) essentially gives 
Yi oHq for all the transitions involving nuclei and is 
thus the same for all of them. 

(b) Double Resonance - Weakly Coupled Spin Systems ;- Eor 

the double resonance case 60- enters the expression for 

^ .1 

(Eq. 11-62) also through the frequency offset of 

the irradiated nucleus and for a weakly coupled spin system 
Eq. (II-62) becomes: 




= (/\E^)i(r± + 


h 

2n 




a A. 


(11-63) 


Thus the contribution of inhomogeneity to the linewidths 
of different transitions is different and the fractional 
change in this contribution (with reference to single 
resonance) is given by the slope of the curve 
versus A^. Any transition for which this slope becomes 

o 

(-Y^/Y^)f ■feke inhomogeneity contribution becomes zero and 
the only width that remains is due to relaxation and 
instability. This situation has an interesting similarity 
with spin echoes in which the fanning out of the spins 
due to inhomogeneity of the magnetic field is completely 
recovered. This also suggests the possibility that by 
double resonance the contrihution of inhomogeneity to the 
linewidths in high resolution n.m.r. can be estimated and, 
therefor^ an estimate of relaxation widths laay be obtained. A 
quantitative consideration of this effect leads to an 
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accurate estimate of the relaxation parameters "by double 
resonance. 

E stimate by Do uhl e Reso nance;- Eor any two transitions 
a — ^ a' and in the double resonance spectrum of a 

weakly coupled spin system the difference in the linewidths, 
from Eqs. (11-63) and (ll-60)j is given by; 


- 'A^>)pp = - Rpp.pp.) + (YiAH,) 




i. 


. "^i 


\\ 

1 

Y- 

2 % 

i 

■^7 4* 

'^3 

2 % 

\ 


AH^ may be readily estimated from this equation by choosing 
two transitions a and p — ^ p ^ having widely different 

slopes provided the relaxation widths are small compared 
to inhomogeneity contributions or that the tvro transitions 
a oi‘ and p — > p ^ have nearly equal relaxation widths, in 
which case the first term in Eq. (11-64) v;ill be small 
compared to the second, and may neglected. The most 
convenient method is to look for, in the double resonance 
spectra, transitions having slopes opposite in sign and 
magnitudes close to unity, formally in double resonance 
spectra transitions v\rith slopes opposite in sign are present. 
As shown below, in weakly coupled spin systems, there are 
always pairs of transitions whose slopes are equal in 
magnitude and opposite in sign. 

Consider a weakly coupled spin system of the type 
Aj^X^. The diagonalized hamiltonian in the rotating frame 
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for irradiation of X spins can "be written as [28]: 




WO' 


(11-65) 


n 


m , 


where I (A) = ^ I„(i) and I„(X:) = which 

i=l . 3=1 . 

i and j are respectively A and X types of spins. J is 

the spin-spin coupling constant between the A and the X 

nuclei. Transitions belonging to A spins are those between 

states |ni(i), m( j)) jm(i) + 1, mC])} . Consider a pair 

of A transitions between the states |m(i), m(3)^> |m(i)-l, 

2i(d)/ and 


m 


(i), -m(3)^ — ^ |m(i) -1, -111(3)^ . Their 


transition frequencies are respectively ; 


. J =>(1))^ . v|^]* m(3) 

-[%A-^ + J (m(i)-l)}'^ + m(3) ^ 


( 11 - 66 ) 


and 




= 27 i;^A_^-[(A^ + J m(i))^ + 111(3) + 

+ [ ^; A-y + J((ia(i)-1) "f- ^ + v^y]"^ m(3)\ . (11-67) 


The slopes of these transitions are given by; 


1 (A^+Jm(i))m(3) [A^+J(iii(i)-1) ]m( 3) 

^ 3^ [(Aj.+Jm(i))UT|j]* " [{Ajj+JCm(l)-l)y“ + 

( 11 - 68 ) 
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and 

I ^ p'’ = _ Ji. 

271 A-^ 2tz ^ A^ 


(11-69) 


and are equal in magnitude and opposite in sign. It is 
then possible to estimate the inhomogeneity contribution 
to the linewidths from the values of these slopes and the 
observed difference in linewidths. 


lie SOLUTIONS FOR X 

11 If 

1. High Amplitude of Irradiation; Bloch Approximation 

Bloch proposed .. [4] that if, in the rotating 

frame in a representation in which the stationary part of 
the spin hamiltonian is diagonal, the difference between 
any two energy levels is large compared to the corresponding 
relaxation matrix elements such that : 



( 11 - 70 ) 


for each a ^ g'^, then the off-diagonal matrix elements of the 
stationary part of the spin density-matrix ( X + *^ 0 ^ 
negligible in that representation. This condition (II-70) 
was shown to obtain even when the irradiation frequency 
coincided with a single resonance transition [9], provided 




(11-71) 


Irradiation strengths termed as high* in this study refer 
to amplitudes of irradiation which satisfy the condition 
(11-71) . 
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✓V 

Under the approximation that ( X + is diagonal 
in the a,a^. ... hasis. (known as Bloch approximation), the 
diagonal elements of 'X are given hy (from Eq. (II~42)): 


2 . 







( 11 - 72 ) 


and the off-diagonal elements of X equal those of (-cr^) . 

Ih^ solution for X obtained from Eq, ( 11-72) depends only 
on the relative values of various relaxation matrix elements 
and not on their absolute values. Thus the double resonance 
spectra for which the condition (11-71) is satisfied, give 
information only about the relative values of relaxation 
matrix elements. 


t% f| 

2. low Amplitudes of Irradiation 

If the amplitude of irradiation does not satisfy 
the condition (11-71) , the entire set of Eqs, (11-40) to 
(11-42) need be solved, in general. These solutions then 
depend on the absolute values of the relaxation matrix 
elements. Two special cases may further be considered. 

If the strength of irradiation is smaller then the linewidths 
in the single resonance spectrum, only intensity changes and 
no splittings are observed in the double resonance spectrum 
(generalised Overhauser effect) [4]. The strength of 
irradiation is termed as "low"^ for such cases and it is 
convenient to obtain a steady-state solution for X('t) 4^ 
the laboratory frame rather thQn in the rotating frame. 
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Secondly, if in tlie double resonance spectrum_ line splittings 
are observed and tlie amplitude of irradiation is not Mgh 
enough to satisfy the inec[uality (11-71) [this is an 
intermediate range of strength of irradiation when 
observed linewidths] , it is possible to obtain a solution 
, in the rotating , without solving the 

complete set of Eqs. (11-41) and (11-42) by mating what 
may be called a partial Bloch approximation [9]. These two 
methods of solution will be discussed briefly in the 
following. 


(a) laboratory Coordinate Solution; - The solution for 
X(t) is obtained by assiming that the diagonal elements of 
X are independent of time and that the off-diagonal elements 
vary as: 



where a, b .... is the basis in which is diagonal. 

The equation of motion of the spin density matrix to be 
solved for this case is given by Eq, (11-14). If in the 
single resonance spectrum there are no overlapping transitions 
present ( line -separations >)> linewidths) ( simple lino 
approximation [rof. 7, p, 522, 15]), it can be assumed 
that the irradiation sith low amplitudes, induces 
only one off-diagoiiol matrix element ^ ^ bctr/e-btL the 
pair^tates affected by the irradiation. If the frequency 
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of irradiation coincides witli a single resonance transition 
a— >13, one obtains, putting T](t) = 0 in Eq. (11-14), for 
the matrix elements of "X('fc) "tiie following equations 
[ref. 7, p. 522, 15]: 

7 \ab'-b' ( ia - Xb' ) = 21 , (II-74A) 

irj a 


2. %ba'a'(^b-^a') , (II-74B) 

a^/b 



(II-74C) 


(11-75) 


where I s Im (zd^) ^ = ^aj ^ja^ , d = 

(^1^2+ ! ^ ^ (bjz'''|a) = ^aj z~ jb^ T 

Equations (11-74) and (11-75) constitute (lT+1) homogeneous 

linear simultaneous equations in E unknowns ( and have 

a unique solution which can be represented as 'Y ;= T /I, 

^'-a a ’ 

where 1^/ are functions of "byp-e of elements and 

are independent of the strength of the irradiating r.f. 
field. The off-diagonal matrix element of X and the strength 
of irradiation are included in I. A solution for % leads 
to^ f or any transition a.' — ^b^j 
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qid^lty^o ^g^/-b^ ^2 

1-^ HI ^2 ' 


(11-76) 


where ^2 “ "^/^hah ^ ~ 

The expression for the intensity of any transition 
a-'" Id can he obtained in a manner similar to that outlined 
in Sec« IIA, hy assuming that the diagonal elements of r|(t)ATe 
zero and that the off-diagonal elements vary as exp (+icJ^t). 
The expression for integrated intensity is then obtained 
as: 


J a'b' 




'o'b' b 


0' a'^a 


(11-77) 


from which the pealc height can be obtained by dividing by 
the total width of the line. 


Por the purpose of the analysis of double resonance 
spectra it is useful to compute the ratio [S/(S’^-S) 

3/ O 

vjhere S° and S are respectively the single and double 

I I / 

resonance integrated intensities of the transition Ou-^D . 
This ratio is given by: 


/-A- 

\s°-s 



[(a 


o^b b 




f 


1 


T 


a'b 


T, 


1 ^ab - ^a^b^ . 


( 11 - 78 ) 
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Equation (11-78) gives a straight line when the ratio 
[S/(S°-S) is plotted against 1/ [ d^ | , having a intercept 

v^'hich depends only on the relative values of various 
relaxation matrix elements and a slope which depends on the 
absolute values. Kuhlmann and Baldeschwieler [15] studied 
the behaviour of this ratio for a range of values of ( d^ | , 
chosen such that line splittings are not observed in the 
double resonance spectra. 

Eor the above plot the absolute change in intensity 
of a transition between two different spectra (single and 
double resonance) is needed. This presents experimental 
difficulties unless there is a reference line in the 
spectrum which is not perturbed by irradiation. To avoid 
thiSj a comparison between the ratios of any two transitions 
a'^ — » b'^and cd^d"^ may be obtained, both in single and double 
resonance spectra. If R represents the ratio of the double 
resonance intensities of any two transitions a^ — ^ b’*' and 
0 .' d'" divided by the ratio of their intensities in single 
resonance: 


R = (s^,^Vs 


c^ d 




(11-79) 


then 


R+1 




a'b -"c^d' 




tVb' - 


(11-80) 
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This equation again gives a straight line having intercept 

^~^^ab d* which depends only 

on the relative values of various relaxation matrix elements 

and a slope 2 /T 2 (T^ ^ -T^ ^ ) which depends on the absolute 
values^ 

If there are overlapping transitions present in 
the single resonance spectrum ' . more than one off-diagonal 

elements of '}C can become simultaneously significant and 
the procedure becomes involved (see Kuhhmann and 
Baldeschwieler [15]}, 

g..artial Bloch Approximation :- If splittings are 
obtained in the double resonance spectrum, the rotating 
coordinate analysis has to be used. However, if the strength 
of irradiation is not large enough to satisfy the inequality 

'transitions ( II~70) will still be true for 
all pairs of states other than the two connected by irra- 
diation, if the other transitions have a separation, from 
the irradiated one, large compared to v^. In such cases 
an approximation that ( % + o'^)has off-diagonal elements 
only between those two states can be made [9], and instead 
of solving all the equations (II-40) to (II-42) , a 
solution for can be obtained by solving only (11+2) 
equation (N equations from Eq. (11-42) and one from each 
Eqs. (11-40 and 11-41). The solutionsfor 'JC will, in this 
case, be dependent on the absolute values of relaxation 
matrix elements. In the limit of no splittings this solution 
reduces to the laboratory frame solution. 
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III - RELAXATION lEBCHAJTISMS MD GOERELAIIOITS 
IIIA INTERACTION HiOOLTONlANS AND CORRELATION CONSTANTS 

TNe liamiltonians for various interactions wMcii can 
contribute to relaxation of spin ^ nuclei in liquid samples 
and expressions for correlation constants are given in this 
section, A general expression for the relaxation hamiltonian 
considered as a random function of time, has already been 
given in Eq. (II-ll) and some properties of spin operators 
A^*^^ and lattice functions Eqs, (II-*12) and (11-15). 


External Isotropic Random fields 


This mechanism generally includes various inter- 
molecular and intramolecular interactions for which explicit 
calculations are not feasible and a combined effect of all 
such interactions is treated as an isotropic random field 
at the sites of the spins of interest [ 29 ]. The sources 
for these fields are external to the spin system of interest 
but not necessarily external to the molecule. The interaction 
of these fields with a spin ' ±'' can be expressed as [ref. 7, 
p. 508]; 

Y. ^ (III-l) 

i 



such that the spin operators and lattice functions as defined 
in Eq. (Il-ll) can be expressed as; 


(+ 1 ) 1 + 

Aj^” = T — i-'Ci) 5 
1 f2 


a[°^ = l^(i) 


(III-2) 
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+ i H|(t)) = 'YiHj(-fc) 

(III-5) 

These random fields are considered isotropic, which is 
introduced through the relation [ref. 7, p. 508, 9]; 




I =Ii/|H+(t)hN 

av 2\Pu^^''| /£ 


7(i^)| 


av = ; 


(III-4) 


where f(i) is the mean square field at the site of spin "* i^'’ 


The degree of correlation "between random fields 
at the sites of different spins in a molecule may "be 


defined as: 


Cii, 




[<hi (*)| 


>av]* ’ 


(III-5) 


where 0^^/ is the correlation constant "between nuclei i and 
i'^ and varies betvreen 0 and l"for the limiting cases of no 
correlation and complete correlation betvireen the fields 


respectively. This definition of is chosen here such 

that the mean square values of the random fields 


(t) 


and / (t) I' 


need not "be equal for 


partially or completely correlated cases as has been 
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considered Toy earlier workers [ref. 7, p. oOS, 11,14-16], 

Extrene-narrov/ing assiunption is generally used for 
calculations witli tliis rieclianisDi, but depends uj)on tbe 
specific interaction being considered as a randorii field, 

2 • Internal Dipole- Dixc ole Interaction 

Tile interaction hamiltonian for direct dipole-dipole 
coupling between different spins can be written as [ref. 7, 

P. 289, 15, 50] 





(III-8) 


where the subscript IT refers to a pair of spins 1^^ and 

is the distance between them, and y| are the spherical 
harmonics of ranlc 2, v/ith 6^^. and specifying the orienta- 
tion of r^^ relative to the laboratory— fixed coordinate 
system. Transforming pl^s to a molecule fixed coordinate 
systemj the correlation constants between any two pairs U 
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and can be obtained by averaging over all the molecular 
orientations 5 as [15>30]; 





( 


r ^ ^ 





20 


(3 cos^ > 

(III-9) 


where is the angle between r^x and r^jx- . Equation 

(III-9) is obtained by using addition theorem of spherical 

Ak 

harmonics. For II = = 0 and g^^^- are Imovm as auto- 
correlation constants and for h called cross- 

correlation constants which can be calculated from the 


geometry of the molecule. 


3» Scalar Coupling 


The scalar coupling between a spin I = ^ nucleus 
and a spin S ^ 1 nucleus acts as a relaxation mechanism for 
spin I, if the spin states of S are having rapid time 
dependence due to e.g, its fast quadrupole relaxation rate. 
The relaxation hariiltonian for scalar coupling can be 
expressed as [ref. 7, p. 310]: 




I(i) • S(t) . 


(III-IO) 


1 

In this the spin S is treated as part of the lattice. The 
sum over ^ i Is over all the spins I. The spin' operators 
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are expressed as in Eq.. ‘(III-2) the lattice operators 

can be expressed as (see Eq,, II-ll) : 

p" (t) = + i- S^Ct) ; P^°Mt) = S^(t) . (Ill-ll) 

n 

There is a complete correlation between the 
interactions of various spins I(i) with the same spin S(t) 
in a molecule and it may be assumed that there is no 
correlation of these interactions with that to another 
spin S'^(t) in the molecule. 


4 . Spin>-Rotation Interaction 


The interaction of the nuclear magnetic moment 
with the magnetic field produced at its site by the rota- 
tion of the molecule containing the nucleus is known as 
spin-rotation interaction and can cause nuclear magnetic 
relaxation. The hamiltonian for the spin rotation interac- 
tion is given by [31], 



(III-12) 


where J is the angular momentum of the molecule containing 

the spin I and C is the spin-rotation tensor. Assuming a 

spherical model of the molecule containing nucleus I, 

one obtains in the molecule fixed frame = Cjt and 

C = C =0, and the interaction hamiltoniani" in the 
X y -i- 


A general expression for asymmetric rotor is given by 
Huntress [32]. 
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molecule fixed frame as [ 31 ] ^ 


where the cuporscriyt ( i) indicates molecule 

fixed frame 

values, and 



A^°^ = ; A^" ^ = T -i- 

n 

+ 

I" , 

(III-14) 

J ^ + 

-{2 ■" 

i J^) . 

(III-15) 

and 



Co = C|t i C+1 = • 

■u • 

(III-16) 


'Transforming A and J from molecule fixed frame to 

r*Ji 

laboratory coordinate by using the relation: 


(A 


) 


'■’■I*'” 

1 


D 


( 1 ) 

'q.q.' 




(III-17) 


and a similar relation for J, -fcKe interaction hamiltonian 
in the laboratory frame is given by [31]: 


<1 


(III -18) 


where 

p(-l) 


(-1) 




D 




(hT_) J 


(l^) 


1 '' 


(III-19) 
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The spin-rotation interaction is completely corre- 
lated for various spins in the same molecule and completely 
uncorrelated for spins in different molecules [31]. 

The spin-rota.tion constant C is usually smaller 
for protons compared to nuclei like and is usually not 

an effecient relaxation mechanism for protons [33v34]. 


5. Anisotropic Chemical Shift 

Under the influence of an external magnetic field 

H the induced precession of the electrons cause«a inter- 

action of the type -v I * A • H , where A is the chemical 

^ 

shielding tensor of ran]-: 2 [ref. 7, p. 315]. The trace 
of this tensor gives rise to the chemical shift and the 
anisotropic part a ' of this tensor is time dependent in 
liquids “because of molecular tum'bling and causes relaxation 
in nuclear magnetic resonance. The interaction hamiltonian 
is given in ref. 7, p. 315. “ 

Uo examples are laiown in which this interaction 
makes significant contribution to proton relaxation, 
prestxmably because of tlio small chemical shifts of protons. 


IIIB COUDITIOUS UITDER WHICH VARIOUS IHTERACTIOIS CAR BE 
APPEOXIMATEB BT A EAlIBOM EIELB 

The random field mechanism is often convenient to 
use, for^ description of relaxation effects, in cases like 
the analysis of double resonance spectra. This mechanism 
is linear in its spin variables and possesses isotropy, 
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resulting in a considerable simplification in tlie calculation 
of a large n-umber of relaxation matrix elements involved in 
the analyses. This mechanism; serves as a model for 
describing the effects of a complex set of intermole cular and 
int^ramolecular interactions often present in a complicated 
spin system, A loiovrledge of the values of the random fields 
at the sites of various spins in a molecule and the correla- 
tions between these fields could then be used to indicate the 
various interactions responsible for these fields. It^ 
therefore, becomes important to examine the conditions under 
which various interactions can be approximated by a model of 
external isotropic random fields. 


1, General Properties of Random Fields 

The hamiltonian for external random fields and the 
properties of their spin and lattice operators, including 
isotropy (Eq- III-4) and the secular approximation (Eq. 11-13) 
have been given earlier. The processes responsible for time 
dependence in the random field hamiltonian will be assumed to 
be Markoff ian, leading to a correlation function G{'r) as 
[ref. 7, p. 297]: 


G(t) =<jE^‘^^(t) F^^^(t-T) 


av 




A, 


(III-20) 


where Tq is the correlation time and is the characteristic 
time for molecular reorientation, which is responsible for 
time dependence of 
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StilDstituting tliG form of from Eg. (II-ll) , 

for a single spin, into the expression for spectral density 
of the “non-viscous liguid^^ approximation Eg. (11-34), and 
using Egs. (III-2) - (III-4) and the relations; 


-1 jfx ^ T 

e e = I'l e 


-i LJij- qx 


(III-21) 


+1 ' + 

where > I” =1“ and 1° = gives: 


+1 


rUV 

apa^ p" 




2t 


g=-l 


[1+(U>^ gr^)^] 




r 




+1 


(-l)^(^jp('l)(t)|2^ 


av 


1 


la''’ 


(III-22) 


In the extreme narrowing case when ^ 1, (see 

Eg. 11-35): 


+1 


jBh 


= 2fT. 


g=-l 


TppTT 


I-lja') 






(III-23) 
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Equations (III-22) and (III-23) for the spectral densities 
contain the essential features of this mechanism. The 
expressions for T^ and T^ can simply he obtained from these 
equations hy calculating the values for a single spin ^ case 

and noti^ that 1/T^ = 2 1^212 ^/'^2 ~ '^1122 ''' '^1212’ 

where I^jl^ = ■§'|l^ and 2^ = -|-j 2^ , 


2 • Var i ous Interac tions 

The procedure for obtaining the conditions under 
which Various interactions can bo approximated by an iso- 
tropic random field will be to obtain the spectral densities 
under '’non- viscous liquid and 'extreme-narrowing ^ approxi- 
mations and compare them with Sqs. (III-22 and 23). All the 
interactions mentioned in Sec. IIIA will bo considered except 
anisotropic chemical shift, which obviously cannot be approxi- 
mated by an isotropic field. 


(q-) D ipole-Dipole Interacti on;- This is a direct magnetic 
interaction between two spins and depends both on the distance 
between the spins and the orientation of intemuclear vector 
with respect to a space fixed frame. The two spins may either 
belong to the same molecule((i) intramoleculan dipole- 
dipole interaction) or different molecules( (ii ) inter- 
molecular dipole-dipole interaction). The time dependence 
in these dipolar interactions arc due to molecular rotation 
in (i) and due to both translation and rotation in (ii). We 
will consider both these interactions separately. 



(i) Intramolecular dipole-dipole interaction ;- If lotli the 
spins* of interest are spin -g- nuclei then this interaction 
is generally treated explicitly in rel$.xation calculations 
through what are Imoim as internal dipole-dipole interactions 
the experessions for \ 7 hicl 1 have heen given in Soc, IIIA, and 
these interactions cannot he approximated as random fields. 
However, if one of the spins has rapid time dependence in its 
spins states arising from a fast quadrupole relaxation rate 
(spin 1) then this interaction, as shoTrti heloTj,can he 
approximated as an external random field. 

The dip ole -dip ole interaction hai)iiltonian Eq. (III- 
6-8) can he rev/ritteii for a single pair as, 

^ (t) = I‘^[S'^ i S" + I^[S"‘'E^-^^ 

+^1 - S" ] 

4“1 +1 

jQl \ qQ. gP (ill — 24 ) 

q=-l p=-l 



whe re 


^“^I /^“*I ^ r<**I — n ^ 

+1 ~ ^-1 “ ^0 “ ^+1 ~ ^O ~ "^-1 

n+l - n”! - 1 /^ . pO _ /s" 

^-1 “ ^+1 ~ “ 4\/3 » o ~ sj3 


1 


(III-24A) 


Substituting this hamiltonian in the expression 
for spectral densities (Eq, 11-34) for *^non-viscous liquid'*' 
approximation and mailing use of the relations (III-21) and 
similar relations for S, one gets; 
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(III-25) 


Consider now tlic case vdien one of tlio spins, say 
spin has a rapid time dependence arising from interactions 
which do not involve I. This can arise from rapid quadrupolar 
relaxation of spin S, if 3^1. The motion of spin S may 
then be assumed to be un.affected by its interaction v\ritli spin 
I and the rapid relaxation of S may be assumed to be 

described by Bloch equations. This condition is usually 
obtained for S = 1 and for S>1 if "‘extreme narr ovring is 
valid for molecular motioiic responsible for relaxation of S 
[ 3]. This latter condition also leads to the 

characteristic times in Bloch equations. In addition if 
3-nd Tp are very much different from t , the correlation time 
for molecular rotation, it is permissible to separately 
a,vorage the Pjand S^in Eq. (I II -25) giving for s ['/]■■ 


Spo) 


s 


2 


w) 


av 


S(S+1) 

3 


HtIA^ 

e 


/s 




(o) S (v 


4” 


0 


av 


= 2S(S+1) 


+itk>gT 


(III-26) 


Equation (III-25) will now have a correlation constant given 
by: 
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1 


X 


/ 

c 


_ 1 _ ^ 1 
'^c '^(}PI+1) 


(III-27) 


In the limit x^ ]_) » true for most of the cases, 

S' Llaking use of this and of Bqs. (111-20,26,27) 

and (11-13), and substituting the value of from Eq, (III- 

P 

24 A) in Bq. (III-25) one obtains for spectral densities: 


^nv 


/ 


2S(S+1) 



l+( CO^qr 


c 



fjipTJ 


X 


< 





(III-28) 


which is identical to the expression (III-22) T;ith 


j _ 2S(S+1) 

^ \~? 

Ihus the conditions for intramolecular dipole-dipole 
interaction to be approximated as random field are that, the 
time dependences of siDin S and that of B are widely 

r>/ 

different and that Bloch equations can be assumed for motion 
of spin S. 


(III-29) 


Correlations :- If there are two or more spins I(i) in the 



molecule having dipolar interaction with a spin S, then the 
correlation between the random fields produced by the 

^ spin at the different spins ^(i) are given, by, using 
Eq. (III-9) and Eq. (III-5); 



(III-30) 




= -i (3 - 1) , 

where '\|/'jjtj/ is the angle between the vectors joining i and i^ 
snins with the snin S, In such a case it is then possihle 

nf 

that the randon fields nay he completely correlated (if the 
three spins lie in a straight line) having different mean 
square fields, Ihe correlation with more than one spin 
having fast time dependence will a.lso he given by Eq. (Ill- 
30), since t^« 

(ii) Intermolecular dipole-dipole interaction ;- In this 
case both rotationa.1 and translational motion of the molecule 
are responsible for changes i3i orientation and magnitude of 
internuclear axis, Hov/ever, assuming in the follov/ing that 
the molecules are spherical and that the spins of interest 
are at centers of the molecules (in vdiich case the rotational 
motion will not be effective ) , expressions for spectral 
densities will be obtained. Ihe case in which the spin is 
not at the center of the spherical molecule does not alter 
the arguments of this section, and will be briefly discussed, 

The interaction iiamiltonian can be expressed as 
( see Eq, III-24) : 


I I b" 


1 


P 


(III-31) 


vdiere the sum is over all j nuclei external to the molecule. 
The spectral densities Eq. (III-25) in this case will then 


become ; 
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where N is the nmlier of spins (or molecules, for this case) 
per unit volume, is the radius of the assumed spherical 

molecule and D is the translation diffusion coefficient. 

The translation correlation time is then given hy x^ = 
2a^/D. Equation (III-33) again represents a random field 
(compare Eq. III-25) having a mean square value; 

f = n|±li S" (III-34) 

0 1 b 

Thus under ext r erne -no.rr owing approximation, and 

if the spin sta.tes of S_ con he assumed to he equally popula- 

z 

ted (high tempera.ture approximation), the intermole cular 
dipole-dipole interaction can he a.pproximated hy a random 
field. 

This conclusion remains unaltered if the spins are 
considered to he situated at a distance h from the centers 
of the spherical molecules. In this case the average distance 
of closest approach between the spins becomes different from 
the above value (2a) and further both translational and 
rotational diffusion contribute to the time dependence of 
this interaction. These two are compensatory effects. The 
correlation function calculated for this model, for b ^ a/2, 
by Hubbard [35] leads to an expression similar to (III-33) 
with a slightly modifie'd value of mean square fields which 
are given by; 

f = liS+li ^ ^2 ^2 ■|^2 ^_^q^233(|)^ + 0.15(|)'^ + ..o..] 

5a’^ 


(III-35) 
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Tile second and third terns add to less than 7 per cent of 
the first ‘term for h = a/2. 

Correlations ; - If there ts more than one spin in the 
molecule, the interaction hamiltonian (III- 31 ) will hecome; 

Z X (III-36) 

i~ d' 1 ?“ 


where i stands for all nuclei in the molecule. The various 
terms -tKoir will occur in the calculation of spectral 
densities will he, in general, of the type; 










V(p+<i)(t) 


(III- 37 ) 


The terms for T;hich 1 = 1 ^^ and d = represent the 
completely correlated terns, included in the above calcula- 
tions. The terms for which j / represent cross-terms 
between dipolar iiiteractions with external spins and may 
be assumed to be completely uncorrelated and taken as zero 
[35] (these terms are zero even when 1 = 1^^, but 3/3^ 
and it is this part v/hich was assumed zero in Eq. (III-32)). 
The terms for which i / i'^ but 3 = 3^ represent the field 
produced at the two spins i and i'' in the i^olecule by the 
same external spin 3 and represent the case of partial 
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correlation. IToggle [14] lias performed a calculation of 
this term hy assuming only transla-tional diffusion and 
neglecting the effects arising spins oeing not at 

the centers of assumed spherical molecules, which also leads 
to f(i) = f(i'^) * It is further assumed that the surface of 
closest approach is anellipsoid of revolution about the 
intemuclear axis of the tv;o spins i and i’^ having a semi- 
major axis (r + h) and semiminor axis h, where 2r is the 
distance between the spins i and i'^ and h is the nominal 
distance of closest approach of spin j to any of the spins 
i (assumed equal for both). The calculation showed that for 
h r, 0 = 1 and for r h, 0 -^0. 


(b) Scalar Coupling :- Using the hamiltonian (III-IO) and 
Eqs. (111-2,11) for the scalar coupling of a single spin I 
Vifith a spin S having fast time dependence in its spin 

fV/ 

states due to fast quadrupole relaxation, and using Eq, 
(III-21), the spectral densities for this case are obtained 


as: 


nv 


400 


471;^ J 



-oo 



iCJj qx 

e _____ 


( |<i(+i)(|a'l+i) 


"i(t) 


/ ^ 

3“^^ (t-x) 


av 



(III-38) 


Assuming Bloch equations for the time dependence of spin S 
(ref. 7, p, 310 and [3]), one obtains q = -q'' and : 



67 


.nv _ 4 tc^J^S(S+ 1) Y' Iql+l) 1 

' 3 V 141+1)3 





(III--39) 


Using extreme narrowing (oJqT^)® 1 where is the 

correlation time for molecular reorientation which causes 

the fluctuations in the electric field gradients at the site 

of spin S (this assumption has to he made for Bloch equations 

to hold for motion of spin S for S')!, [3])^ we obtain = 

= T where r is the quadrupole relaxation time of S, 

^ Q. Q. 

Equation (II 1-39) then represents a random field (compare 
Eq. 11-22 and note (hJ^ instead of in the denomi- 

nator)witha mean square value: 


f = (III-40) 

Thus the only condition for scalar coupling to he 
approximated as a random field is extreme narrov/ing f or 
molecular reorientation, which is generally true. It may, 
however, he noted that (CJ-^- - "''q normally not 

less than unity and in fact quite often ( - OJq)^ ^ 1 
is obtained. 


Correlations; - If there are number of I(i) spins coupled 
to one S spin (having fast time dependence in its spin 


states) then the random fields at various I(i) spins will 
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iDe completely correlated (C = 1) witii mean square fields 
equal to „S(S+l)/3. 

In case there are more than one spin of the type 
3(3) (having fast time dependence) the portions of the 
fields produced at the site of any spin I(i) hy different 
SCj) spins may be assumed to be completely uncorrelated 
v/ith each other. This then leads to a case of partial 
correlation between the fields at various I(i) spins with 
the extent of correlation governed by the coupling constants 


(c) Spin-Rotation Interaction ;- Using the hamiltonian 
(III-18) and Eq, (III-21) the spectral densities are 


obtained a si 
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This expression has been evaluated by Hubbard [31], under 
the assumption that the characteristic time for the 
change of the components of angular momentum is very 
much different from ' ' 2 molecular reorientation, as: 


1 

TJqpry 


ikT y 

9'fe^ 


2(C„-Cj.)^2Tj^g 

[l+(Wj 



(III -42) 
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v/liere 



(III-43) 


and either '^2 n-umher 

of identical I spins on the same molecule and hi is 
the Boltzmann factor. In y/riting expression (III-42), 
no correlation hetyreen spins on different molecules and 
complete correlation between spins on the same molecule 
(Mre assumed [33]. 


Eq, (III-42) represents a sum of two random fields 
having different correlation times and different mean square 
fields. However if 

a random field with mean square field as: 


<< 


<'^2 
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^ ^1* ( 111 - 42 ) gives 
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f = If 
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(20,,^%)' + 2(Cjj - Cj_)^ 


= 1 


Ihl 


(2Cf+C|.) 

(III-44) 


I lie TEEAlMEHl OE SCAL.1R COUPLING AS A EELAXAlIOlf LIECHAlfISM 


When a spin -g- nucleus is coupled hy scalar coupling 
to a nucleus having spin 1, this coupling generally does 
not appear as a splitting due to the rapid quadrupole 
relaxation rate of the latter. Depending on the relaxation 
time T of the quadrupole nucleus and the magnitude of this 

Si 

coupling a broadening in the spectrum of spin ^ nucleus is 
generally observed. When this coupling is entirely * washed 
out"^ [ 37 ] (no residual splitting remains and only broadening 
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is observed) the scalar coupling interaction is treated as a 
relaxation mechanism. In the following the conditions under 
which such a treatment of scalar coupling becomes valid are 
(M-scussed. 

The manner in which the coupling of I (I = -I) with S 
(S affects the single resonance spectrum of spin ^ can 

be treated by considering the total hamiltonian: 
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(III-45) 


whe re 


and ^^0 zeeman terms of spins I and £ 


respectively, scalar coupling between them, 

is the interaction of the quadrupole moment of ^ 

with fluctuating electric field gradients and is responsible 

for relaxation of spin S, "^(t) is the interaction of 

■' /)/> 

(t) 



observing r.f. field with the spin system and ^ 

represents mechanisms of relaxation of I other than scalar 

coupling with S. The equation of motion of the spin density- 

* 

matrix a can then be written asj 


do 

dt 




o 



sc 



(t), <l] 




( 111 - 46 ) 

This equation includes complete information regarding the 
dynamical behaviour of the spin system including spin S. In 
order to describe the relaxation effects on a transition 
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(m(i) in. the spectrtmi of spin I, one considers 

simultaneous existance of off* diagonal elements of a of the 
type |m(i), m^^ 4 — > jm(i^), 1°^ a-11 different values of 

m , the eigenvalue; of S . furthermore, as may he seen from 
Eq. (111-39) that if 47 r^ - ^ 3 )^ if the relaxa- 

tion time T of spin S due to '^^„(t) is such that 
(uJ^ - yy 1 , the terms like (l.S + I S,) in I*S 

will not contribute significantly to various relaxation 
matrix elements. In such a case oyd.y fhe terra 
retained in With such an approximation and setting: 


, With such an approximation and setting 


^i(t) = 0 in Eq. (III-46) an expression for I = and 
S = 1, for the line shape of the spectmm of spin I, is 
obtained, with the relaxation time of spin S as a parameter, 
as [ 38 ]: 


I(x) 


(45+r,^+5x^p® ) 

2 2 5x^ +r]^ ( 34x^ - 2x^^ +4 ) +ri'^ ( x^- 2x'^+x' 


(III-47) 


where p = IOtct^^J and x = AUJ/amJ, in which is the 

difference between the 1 armor frequency of spin I and the 
observing frequency. Equation (III-47)> given by Ibragam 
[ 7 , p. 504 ], was also obtained through basically similar but 
slightly different arguments by Pople [39] following a 
general treatment of lineshape by Sack [ 40 ] and Anderson [41] . 
This equation describes the spectrum of spin I for the entire 
range of values of If 1, Eq. (III-47) leads to: 
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a lorentzian having a half-width at half-intensity equal to 

J^/3. 

q 

This is also the linev/idth and lineshape obtained 
when scalar coupling is treated purely as a mechanism of 
relaxation with the condition that {tO^ - yy 1. 

From the expression for spectral densities for scalar 
coupling mechanism Eq. (III-39), (assuming extreme narrowing'' 
for motion of spin one obtains for and T2 spin I: 





S(S+1) 


l+(Cl^j “* g ) 




(III-49) 


a.nd 



4Tt^J^S(S+l) 


/ 

^ l+(C\9j- ^ 3 ) ^"^q 



(III-50) 


From Eq. (III-50) it can be seen that, under the condition 

(Mj - S = 1, the linewidth (I/T2) is 

the same as in (III -48) . 

Thus the conditions under which scalar coupling can 
be treated as a relaxation mechanism are that r\^ 1 and that 

extreme— narrowing holds for motion of spins S. 

Equation (III— 49) } shows that if the above mentioned 
conditions i.e. 4m"j2 « (<^j - COg)^ and (40^ - iUlg)^T:^ 1 
hold, the contribution of scalar coupling to becomes 

negligible and in that case the scalar coupling mechanism 
serves only as a source for spin-spin relaxation. The contact 



73 


of the spin system with the lattice is therefore primarily 
through other interactions in the sample. 


HID GROSS-TEEMS BETTOER VARIOUS IRIERACTIORS 


In a sample various interactions mentioned in the 
previous sections often simultaneously co-exist and it is 
only their relative importance that differs from sample to 
sample. In such a case there can exist cross-terms Between 
different interactions in the calculation of spectral 
densities such as; 
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Besides cross-terms within an interaction. Ihe cross-tems 
within an interaction have already Been discussed under the 
title '^correlations'*' in each case and in many cases found 
significantly important. The cross-terms Between different 
interactions normally vanish Because of either very different 
orders of magnitudes of time dependences involved in their 
hamiltonians ,in which case the two parts can Be separately 
averaged giving zero for the averages, or difference in 
transformation laws for each of the hamiltonians yielding 
zero on averaging [ref. 7, p. 309]. Specifically the 
cross-terms Between scalar coupling and either of the dipole- 
dipole interactions or the random fields vanishes Because of 
different orders of time-dependences involved. 



74 


Chapter IV 
EXPERniBSriAL 

Page 

IVA PEBQUEICY SI7EEP DOUBLE ilESOlTMCE SEI-UP 75 

1. Introduction 75 

2. Basic Principle 76 

3. Details of Exp er incut al Arrangement 79 

(a) lock Channel 80 

(h) Ohservation Chainiel 83 

(c) Irradiation Channel 86 

(d) Pield- frequency Control Panel 86 

(e) Operational Details 86 

4. Performance 91 

IVB SAI;IPLE PREPAMIIOiT 91 



75 


IV - BXPEEIMSiraAl 

IVA FEEQUEHCY SWEEP DOUBLE SESOlUUCE SET-UP 

!• Introduction 

In homonuclear double resonance experiments the 
irradiated and the observed transitions belong to the same 
nuclear species, and- therefore the frequencies of the 
irradiating and the observing r.f. fields differ by an audio- 
frequency. In such experiments the second r.f. field is 
produced by modulating the magnetic field with the appropriate 
audio frequency. Double resonance escperiments can, in 
principle, be performed both by changing the magnetic field 
and keeping the frequency fixed (field sweep operation) or 
by keeping the magnetic field fixed and changing the frequency 
(frequency svreep operation). The frequency sweep experiments 
are convenient for relaxation studies and are used in the 
present work. 

The conventional high resolution n.m.r. spectro- 
meters, like the Variaii IIP-100, use field sweep operation, 

Tliis instrument has been modified to perform frequency sweep 
single and double resonance experiments, by using field- 
freqiiency-locking- control loop [23,25]. Single resonance 
spectra have been obtained without loss of resolution and 
with improved stability, compared to field sweep operation. 

The modification, however, does not disturb the field sweep 
operation. The details of this modification are described in 
this section. 
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2. Basic Principle 

In the first stage of this er:periment the magnetic 
field is held fixed at a particular ralue , The exi^erimental 
sample is addedioifLa small amount of tetramethylsilane 
(HIS) \ 7 hich serves as a source of the locking signal. The 
magnetic field is modulated hy an audiO' frequency f^, Fig. 
(lY-1), and is hrought in resonance with the first low-field 
sideband of HIS. The output of the r.f. receiver is then 
phase detected through a loclr-in-ampllfier Trhich is 
referenced by the same audio oscillator and has amplifiers 
turPod to frequency f^^. Hie phase difference between the 
signal and the reference of the lock-in-amplifier is so 
adjusted that, as one sweeps through the resonance, the d.c. 
output of the lock-in-amplifier is obtained in the shape of an 
error signal. This error signal is fed to the flux-stabilisier 
such that it gives an appropriate correction voltage r/henever 
the field drifts in either direction from its locked value. 

The output of the r.f. receiver is also monitored on a scope, 
after amplification through an amplifier tuned, to frequency 
f^^. A continuous steady display of a signal at f^ indicates 
a locked field. 

The single resonance spectrum is obtained by modu- 
lating the magnetic field by a second audio oscillator and 
sweeping its frequency (f 2 ) by a synchronous motor, Wlien 
fg = % ” 6, where 6 is the frequency difference between the 
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Block diagram of the experimental arrangement 
for a frequency sweep homonnclcar douhle 
resonance sot- up. 
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signal to "be observed and IMS, tli8 sideband of the signal will 
come in resonance with the loched field. The output of the 
receiver which nov/ contains a component at f^ also, is fed to 
another lock-in- amplifier which is referenced hy the oscilla- 
torgenerating fg. This lock-in- amplifier is a hroad-hand 
amplifier. The output of the lock-in-amplifi er is fed to a 
recorder. The sweep oscillator has a linear frequency dial 
so that the different parts of the spectrum are all scanned 
at the same rate and tha,t tho spectra are obtained with a 
linear frequency scale. 

The irradiation of the spectrum is achieved by 
modulating the magnetic field vath a third audio- oscillator 
(f^) such that fj = f^ - 6^, where 6"^ is the frequency 
difference between tho line to be irradiated and TMS. The 
amplitude of this oscillator is controlled by a calibrated 
attenuator . 

5 • Details of the Experim enifcal A rrangement 

The high resolution n.m.r. spectrometer HE-100 
consists of a 12 electromagnet (max. field 23.5 x 10''^ gauss) 
havir^ a current regulated power supply, A flux stabilizer 
V3506 controls the time variations in the magnetic field while 
a field homogeneity control unit smoothers out any spatial 

m 

variations of the field in the region of the sample. Crystal 

controlled r.f. transmitter and receiver Y4311 operating at 

100 MHz are used for proton resonance. Tho homogeneity and 

8 

the stability of the magnetic, field are --»^1 in 10 . The 
resolution of the instrument is 0.5 Hz. A linear sweep 
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miit T4352 giving a sav-^tooth is used for scope presentation 
of tlie signal and a slovj sv^eep unit V3507 whicli operates 
in conjunction with the flu3: stabilizer for slo\7 recording 
of the spectrum in field sweep mode, 

(a) Lock Channel ;- The purpose of this channel is to hold 
the magnetic field fixed at a particular value. The 
2 KHz audio-oscillcitor and lock-in-amplifier available in 
Variant V3521A integrator, have been used for this purpose 
fig. (IV-2) . The output of the audio-oscillator is attenuated 
by a G-R 1450-TB attenuator, which serves as an isolation 
stage, before feeding to the svireep coils. Since the lock 
signal is derived from side-band resonance of TlIS, the 
integrator is operated in'* side-band mode. The receiver 

output at J314 . on V4311j which now contains an audio compo- 

% 

nent at 2 KHz, is fed to the signal channel (referenced 
internally) through an attenuator. The phase difference 
between the reference and the signal can be adjusted vmth a 
front panel control. The output of the lock-in-amplifier, 
which is the error signal^ is taken from the ’’’'recorder output’^ 
of the integrator and fed to the flux stabilizer through a 
switch. This error signal may be fed to pin 5 of J1 on the 
sloYj sweep unit. The amplitude and phase of the error signal 
are adjusted for maximum stability of the lock. The stability 
is monitored by feeding the output of the receiver at J314 
to a scope, after amplification through a GR 123 2A tuned 
amplifier tuned to f^. An amplitude of the order of 100 
millivolts is usually needed for the error signal. 



SI 


Block diagram of the lock-chajmel. Iho 
compononts sho^m mthin tho hrokcn lines are 
arailahle in the Varian integrator V5521. 
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("b) O~bser-yation Cliaxoiel i-- A G-E 1304-B, l)eat-freq.uency 
audio-oscillator, eq.uipped with a linear increment dial 
(range 100 Hz) is used in this channel, Hig. (IT-3) . She 
output of this oscillator is fed to the sweep coils at the 
modulation input jach on linear sxveep unit through a HP 350-D 
attenuator. The loch-in-aniplifier used in this channel is an 
' EHG "^(Electronics Missiles and Communications Inc. Mount 
Yemen, USA) model EJB. Since the input impedence of the 
reference channel of this loclc-in-amplifier is 100 K-il, the 
reference signal from heat-frequency oscillator is fed 
through a impedence matching net work. 

Ihe signal channel of the EBiC lock-in-amplifier 
contains a negative feed hack pentode amplifier which is 
heavily degenerative except at null frequency of a tvrLn-tee 
net v7ork, placed in the negative feed hack circuit. Since 
the frequency is sv;ept in the present experiment, this lock- 
in-amplifier is to he operated over a range of frequencies. 
Ihese twin-tee ^s were, therefore, eliminated from the circuit 
and the screen grid of the pentode amplifier grounded, and 
as there is no degenerative feed hack at any frequency the 
circuit acts as a hroad hand audio-amplifier, 

She output of this lock-in-amplifier is obtained 
through a low frequency hand pass filter and fed directly 
to a G-14 recorder. The phases of the frequency sweep 
spectra are adjusted hy adjusting the reference phase of the 
EJIG lock-in-amplifier. 
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Fig. IV- 3. 


Block diagram of the ohserv'ation. chamcl. 
in EMC, model lock-in- amplifier and a 
GR 1304~B he at- frequency audio oscillator 
are used in this channel. 



lOOK lOOK 


8 



OBSERVATION CHANNEL 
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(c) Irradiation Ohairael :- A HP 205-AG audio-oscillator has 
been used for irradiation. Phis oscillator has a maximum 
output pov/er of 5 watts and a calibrated attenuator in its 
output, Phe output of this oscillator is fed at the 
modulation input jacl: on linear sweep unit. 

A HP 5210 frequency counter is used for measuring 
frequencies of various audio oscillators. 

(<3.) Pield-frequency Control Panel: - A control panel 
essentially consisting of various switches, for a quick change 
from field-sYreep operation to frequency sweep operation and 
vice-versa has been installed just below the oscilloscope. 

One switch, Pig. (IV-4), changes the input to the oscilloscope 
Y between integrator 'scope output'^ and the tuned amplifier 
(2 KHz) output for monitoring the lock. Another sv/itch 
changes the recorder from integrator 'recorder-output"* 
to 3L0.. " reccrrder-output ' . One oh/off switch 
is used for the error signal to the flux stabilizer. Phe 
amplitude of the error signal can be adjusted by the poten- 
tiometer Pune lock'^'. 

If all the switches in this panel are flipped 
dovm the spectrometer is in the field-sweep mode, 

(e) Onerational Details ;- Phe magnetic field is adjusted 
for best homogeneity and resolution in the field s\7eep mode 
and the low field side band of PMS obtained in the shape of 
an error signal in this mode. At the mid point of this 
signal the sweep is taken off and the error signal fed to the 
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Fig. IV-4. Diagrajn sliomng the functions of various 
switches Qn the field-freq.uency control 
panel. 
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flux stalilizer. Tlie sign of the error signal is chosen "by 
trial and error ly noting the locking or antilocking of 
the field. The lock is made stable by adjusting the ampli- 
tude of modulation and the amplitude of the error signal 
being fed to the flux stabilizer, A slight adjustment of 
the r.f, reference phase at V4311? bas also been found 
useful for stabilizing the lock. The stability of the lock 
can be tested by deliberately trying to drift the field 
slov/ly and noting that the error signal compensates for it. 
Small variations in homogeneity are adjusted by maximising 
the output of f^ at J314, ?/hich is monitored on the scope. 

kor the purpose of scanning the spectra the linear 
increment dial of the beat-frequency . oscillator has been 
coupled to a synchronous motor through gears. Different 
driving speeds are obtained by changing the gear ratios. 

The amplitude of the be at -frequency oscillator is 
usually kept large for the following reasons. As the input 

impedence of the reference channel of the EMC lock-in- 

ityC 

amplifier is 100 Ki^-, the output offbeat-frequency oscillator 
is first dropped across a 1 K-^'^ resistance to match the 
output impedence of the oscillator and then dropped across 
a 100 Knu resistance to match the input impedence of the 
reference channel of MG (Eig, IV~3) . The two resistances 
are isolated by a 100 Kai potentiometer. The modulation 
amplitude of f 2 and the signal to the reference are adjusted 
respectively by the HP 350-D attenuator and the 100 K li. 
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potentiometer. A large amplitude of the heat~freq.uency 
oscillator then ensures a high attenuator setting and a 
high -value on 100 If-iU. potentiometer. A high value of 100 KiX 
potentiometer avoids any mismatching of the reference input 
imp e dene e . 

The frequency of irradiation is made to exactly 
coincide -with a singl.e resonance transition in the folio-wing 
manner. In the absence of irradiation the observing 
frequency f^, is swept to the peak of the line to be irra- 
diated and the motor is then s-witched off. The amplitude of 
irradiation is then made fairly low (of the order of the 
amplitude of the obser-ving fiSeJ.d. ) and its frequency 
adjusted as near to fp as is possible with a frequency 
counter (within fraction of a Hz) , At this stage the output 
of the EI/lC lock- in-amplifier will beat at the difference 
frequency of these two oscillators. This beat can be 
monitored either at output meter of the ElfC lock-in-amplifier 
or on the recorder. The frequency of irradiation is now 
carefully adjusted by a fine control to reduce the beat 
frequency to zero, at which the IMG lock-in-amplifier -will 
show a large d.c. amplitude as indicated by a large deflec- 
tion of the output meter. 

The irradiation oscillator is normally operated at 
full power and various amplitudes of irradiation^^ol>tained 
by using the calibrated attenuator. This facilitates (if 
all the other parameter are kept constant) a calibration of 
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complete range of strengths of irradiation once one of the 
values is calibrated. 

4 * Performance 

A comparison between the field-sweep single 
resonance proton spectrum of 2-bromothiazole and its fre(i- 
uency— sweep spectrum obtained with the above set-up is shown 
in Pig. (IV-5) . Ihe resolution in the frequency sweep 
operation is of the same order as the field sweep operation, 

IVB SA3,!PLE PREPARATION 

Commercial samples have been used in this work. 

The samples were purified, by recrys tali zing in case of 
solid samples and redistillation in case of liquid samples, 
before sealing in 5 m.m. outer diameter n.m.r. sample tubes. 
A few drops of tetramethylsilane was added in each sample 
to provide the locking signal. The samples v/ere degassed 
under vacuum by a series of freeze-pump-thaw cycles [14,22] 
and finally sealed in vacuum. 



lig. IV- 5. Ihe single resonance proton spec trim of 

2-1 romothi azole (a) in the field sweep 

operation and (h) in the frequency sweep 

operation obtained with the present set-up. 

^ f 

The broader doublet is assigned to A proton 
(position 4). fhc cheinical shift betweegi 
the two groups is i *1^ - = 18.8 + 0,1 

Hz and the coupling constant 
0.05 Hz, 


= 0.55 + 
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Chapter V 

DOUBLE BESOMCE STUDY OF PROTON BELAXAgIQK II A 
SYMMETRICAL THREE SPIH SYSTEM (ABg) 
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DQTJBLB HES0MG5 S'::UDY 01^ P50T01 BELAXATIQH IIT A 
SYiaiETEICAi ■TIIiffiE SPIE SYSSEIi (A3p) 


VA IlISRODUCTION 


In this chapter the results of a proton-proton 
double resonance study of spin relaxation in a syimnetrical 
three spin system (AB 2 ) formed by the ring protons of 2,6- 
dibromo aniline are presented. The sample consisted of a 
20 per cent by weight solution of 2,6-dibromoaniline in 
CCl^. Double resonance spectra exhibit significant relaxa- 
tion effects over a wide range of strengths of irradiation 
(from 0,1 Hz to 9.0 Hz). Ihe effect of magnetic field 
inhomogeneity on the linewidths of double resonance transi- 


tions has been explicitly considered in this analysis and 
is found to be of considerable importance in deriving infor- 
mation on the relaxation processes. Sv/o relaxation mecha- 
nisms ,(i) external isotropic random fields and (ii) internal 
dipole-dipole interactions between the ring protons along 
with possible correlations between them, have been considered 

^ if 

for detailed calculations. The calculations for high 

w u 

strengths and low strengths of irradiation have been 


performed by the methods outlined in Chapter II. A detailed 
comparison of the calculated a for both high and ' low 


strengths of irradiation with experiment is found 


to be essential for distinguishing between the contributions 
of the different mechanisms in this study. 



96 


VB SINGLE RESONANCE SPECTRUli 

RjS-DilDromoarLiline "belongs to a symmetry point 
group 02 ^- and its ring protons form an AB 2 typo of spin 
system (in the notation of Pople et.al. [2]). Ihe resonance 
of proton^' A (in the para position with respect to NHg 
group) is expected to appear at a higher field than that of 
protons ' (meta to NH 2 group) This gives "> where 
h' is resonance frequency. The "basic-functions , diagonal and 
off-diagonal matrix elements of the spin hamiltonian and the 
transition frequencies and intensities for an AB 2 spin system 
are given in Table (7-1), The diagonal and off-diagonal matrix 
elements of double resonance hamiltonian in the rotating 
coordinate are also indicated in this table. In single 
resonance there are six symmetric and two antisymmetric spin 
functions, whose frequencies depend on two parameters viz. 

chemical shift between the a ' and B protons. 

The single resonance spectrum of the molecule is 
shown in Pig. (V-1) . The spectrum has four well resolved ’A 
transitions and two closely spaced doublets in '' B* transitions. 
The small separation in the two ’"'b'"'’ doublets has been 
determined in the field sweep operation by wiggle-beat ^ 

NH 2 is an electron donating group and produces greater 
chemical shielding at the site of the para proton than that 
of the meta proton in mono substituted benzenes, while 
bromine is inactive in this respect [ref. 2, p. 259]. 


AB‘ 


the spin- spin 
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Fig. V-1. Frequency sweep single rcsonanoo spectrum 
of tlie ring protons of 2,6-dit)romoanilinc , 
which form an "typo of spin system. The 
chemical shift between the two groups is 
I V. “ "^-D I = 89.95 Hz and the coupling 

f A I 

constant is = 7.95 Hz. 



100 




102 


Jig. V-2. Energy level diagram of an spin system 

for which ^ spin functions noted 

against different energy levels are the 
respective spin product functions. The 
various transitions are also shovm. schemati- 
cally with heights proportional to the single 
resonance intensities. The lines 1 to 4 are 
transitions and the lines 5 to 8 are 
B transitions. Line 3 is the transition 
between antisymmetric states 7 and 8. 
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linc-widtli changes. She values of Vg for spectra vjhich 
show splittings were determinGd in the following manner. 

■The hamiltonian in the rotating coordinate ^ (Eq., 11-26) 

is diagonalized by a computer for various values of Vg for 
any particulo^r line of irradiation and the calculated line 
positions in the double resonance spectra plotted against 
Vg. These plots are then compared with the observed spectra 
for irradiation of that line [11]. This leads to a deter- 
mination of Vg to an accuracy of<^+ 0.05 Hz, For cases of 
weak irradiation, low values of Vg were produced in the 
experiment by using a calibrated attenuator, and therefore 
for these spectra the values of V 2 were calculated by 
multiplying the values determined for strong irradiation 
cases, by the respective attenuation factors. 

The spectra for low ^amplitudes of irradiation 
exhibit the following general features. Figure (V-3(a)) 
shovfs the spectrum of ^ group while line 6 of group 
is irradiated by an amplitude of 0.1 Hz. This produces 
a decrease in the intensity of line 1 and an increase in 
that of line 2. The corresponding spectrum for irradiation 
of line 7 (Fig, ‘V’-3(b)) shov;s a reverse effect in the 
intensities of lines 1 and 2. Both these features can be 
qualitatively understood by referring to the energy level 
diagram (Fig, 7-2) and noting that the irradiation of a 
pair of levels causes an increase in the population of 
higher level and a depletion in that of the lower level. 
Therefore for irradiation of line 6 (6 4 transition) the 
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lig. V-3. Disserved and calcula1:ed doulslG resonance 

li i 

speotra of A group in 2, 6-di’bromoaniline, 
while (a) line 6 and (h) line 7 of B group 
irradiated with an amplitude of 0.1 Hz, The 
calculated spectra are for random field 
mechanism (mechanism (i)) mth f(A)/f(B) = 
2.0, Og 2 = 1*0, = 0.0 and 2f(B)T^ = 0.2 

rad,/ sec. 
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intensity of line 1 (6 5 transition) should decrease 

while that of line 2 (4 3 transition) should increase, 

in conformity with the observation (Fig, V-3(a)). These 

tiKe 

changes are reflected in^intensities of all those transi- 
tions in the spectrum which have an energy level in common 
A'j'ith the pair being irradiated"^. Similar intensity changes 
have also been observed for weak irradiation of other lines 
except line 5. Typical cases of spectra when lines in 
A group are irradiated vath Vg = 0.1 Hz are shown in 
Big. (V-4). 

A weak irradiation of line 3, which is the transi- 
tion between the antisymmetric states? does not produce 
intensity changes in other lines'^ and furthermore in any 
of the observed double resonance spectra (Figs. V-3,5 and 6) 


f Intensity changes, smaller than the above, can also occur 
in the other transitions \jhich are connected by relaxation 
processes to the irradiated pair. These changes are not 
obvious from an inspection of the spectrim but may become 
significant for understanding the spectmm' on a quantita- 
tive basis. 

Intensity changes and splittings occur for irradiation 
of this line for large values of due to the fact 
that the symmetric states belong to the neighbouring 
transitions (particularly line 2 (4 -V 3 transition) , 
get perturbed by the irradiation. 



108 


lig, V-4-. Observed and calculated double resoiianco 

spectra of B* group in 2, b-dibromoanilino, 

7/iiile (a) line 1 and (b) li'ie 4 of ' a * group 

are irradiated with an amplitude of 0,1 Hz, 

Ihe calculated spectra are for random field 

meohaiiism v/ith f(A)/f(B) = 2.0, 0^. ^ = 1.0, 

J3^jd2 

- 0,0 and 2f(B)T^ = 0,2 rad, /sec. 
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V"5. Observed and calciilated double resonance 

spectra of *A group in 2, 6-dibromoaniline, 
while (a) line 6 and (b) line 7 of group 
are irradiated ?/it}i an njnplitude of 2.2 Hz. 
liio calculated spectra are for random field 
meclianism v;ith f(A)/f(B) = 2.0, C-p = 1.0 
and 0^ = 0.0. 



Ill 
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Pig. V-6. Observed ajid calcLil ate d double resonaiioe 

'.'pectra of "'A* group in 2,6-dibromoaniline, 

while (a) line 6 and (b) line 7 of ’ B* group 

are irradiated vatli an axaplitude of 6.75 Hz. 

The calculated spectra are for random field 

mechanism with i(A)/f(B) = 2.0} Gt, -n — 1.0 

12 


and = 0.0. 
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Fig. V-7. Observed ab-d calculated double resonance 

4 * 

spectra of B group in 2, b-dibromoaniline , 
while (a) line 1 and (b) line 4 of A group 
are irradiated by amplutudos of 5.25 Hz and 
9.25 Hz respectively. Ihe calculated spectra 
are for random field mechanism with f(A)/f(B) 
= 2.0, Cg 3 = 1.0 and = 0.0. 
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tile intensity of this line is not changed^. This indicates 
that the relaxation hamiltonian does not connect the symmetric 
and antisymmetric states. 

2. linewidths 

The double resonance spectra also exhibit significant 
linewidth variations, For example the spectrum of 'a* proton 
obtained when line 6 of the spectrum of B proton is irra- 
diated by an amplitude of 2.2 Hz (Fig, V-5(a)), has two trip- 
lets nearly symmetrical in frequency with respect to center, 
but having considerable difference in their linev/idths. This 
asymmetry gets exactly reversed when line 7 is irradiated by 
the same amplitude (2.2 Hz) (Fig, 7-5(b)). It v/as pointed 
out in Sec. (IIB-2) that the inliomogeneity in the static 
magnetic field can contribute differently to various transi- 
tions in the double resonance spectrum and that this contri- 
bution to any transition a — s a'^ may easily be obtained 
from the plot of 6J //Bm versus A- for weakly coupled spin 

CXCX j 

systems. For making this estimate the present system may be 
considered as a weakly coupled system to a good approximation 
(J/ - Vg! ~ 0.09). A plot of the frequencies of the k* 

transitions - W^)/2m versus A^ = ( OOg ~ /2.% for 

Vg = 2.2 Hz, is shown in Fig. (V-8) . These plots are 
constructed by diagonalizing with the help of the 


f Line 3, therefore serves the useful purpose of being an 
intensity reference in the double resonance spectra. 
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V 8. Plots of oalculatGd doublo resonance frequen- 
cies of *a' transitions versus ‘'frequency- 
offset of irradiation of 'b"’ nuclei. These 
plots are calculated for the spin system 
for v^ = 2.2 H 2 . The numbering of the states, 
corresponding to the transitions, is such 

that for large positive values of and 

s 

Vg = 0,0, those states become the single 
resonance states as indicated in Pig. V-2. 
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computer and de ^jernir-ing the corresponding transition frequen- 
cies for different values of keeping constant. In 
joining the calculated points to construct these plots 
certain non-crossing rules [42] should he oheyedt The tran- 
sition frequencies in a freqiiency- sweep spectrum for any 
value of Ag are given hy the points of intersection of an 
ordinate^ dravm passing through tha.t point, vdtli the contours 
of Fig. (Y-8). According to Eq. (11-63) the change in the 
contribution of the inhomogeneity to the width of a particular 
transition is then determined by the slope of the corres- 
ponding contour at the point of intersection. It may be 
seen from Fig. (V-8) that the frequency of the transition 
between the antisymmetric states 7 8 is not changed_^due 

to a change in A^, giving a zero slope for this line. This 
is a consequence of the fact that the irradiation does not 
mix states of different symmetry. This line (line 3) will 
then have the same contribution of inhomogeneity to the line- 
width as in single resonance. The frequency of the transition 
3 — f 4, which is weak in intensity, changes very little mth 
A-g, since the coupling between A and B nuclei is not very 
strong in this case, In a wealcly coupled spin system AXgj 
this line merges with 7 8 [43] and does not show frequency 


An operational statement of the non-crossing rule is: two 
contours belonging to two transitions which have an energy 
level in common will not cross in the presence of a 
perturbation. 
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shifts. However, in a strongly coupled spin system it 

^ d. 

should show frequency shifts. 

Irradiation of line 6 corresponds to Ag = 3»6 Hz. 
From the slopes of different contours in Fig, (V- 8 ) at this 
value of can be seen that for transitions 5 — 2 , 

5 — > 4 and 5 -^ 6 , forming the lov/-fre quency triplet, the 
inhomogeneity broadening is mostly reduced while it is 
increased for transitions 1 5 ^ 2 , 1 4 and 1 6 , forming 

the high-frequency triplet in the same spectrum, in agreement 
with observation (see Fig. V-5(a)). Through a similar argu- 
ment it is readily seen that this asymmetry is reversed when 
line 7 is irradiated, which corresponds to A-g = -3.9 Hz 
(see Fig. 5(b)). It may further be seen from this plot that 
at these ordinates the slope is nearly^ +1 for various tran- 
sitions for which there is a increase and -1 for those for 
which there is an recovery in the inhomogeneity contributions. 
This gives an inhomogeneity contribution of 2y£\E^ to the 
linewidth of the former and nearly zero for the latter. For 
transitions 3 6 and 3—^2 the slopes at these ordinates 

are almost zero (Fig. V-8) giving the contribution of inhomo- 
geneity to the linevidths of these transitions same as in the 
single resonance. Making use of these observations, and 
the observed linewidths for irradiation of lines 6 and 7 with 
Vg = 2.2 Hz, and the arguments of Sec. IIB-2(c), a value of 
(yAHq)/ix is obtained to be about between O.I 5 Hz - 0.18 Hz 
(full-vddth at half height). 
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5 ' Qaloulation of Relaxation Bffe _ct_s 

llie theoretical double resonance spectra were calcu- 
lated according to the ‘theory sketched in Chapter II for two 
relaxation mechanisms, (i) isotropic external random fields 
and (ii) internal dipole-dipole interactions, the interaction 
hamiltonians and expressions for correlation constants for 
?/hich have been given in Sec, IIIA, The spectra for large 
values of the strength of irradiation ( ^^ 2 ! | aa'^ ^ 

calculated using the Bloch approximation, Bor'low' strengths 
of irradiation, in which only intensity changes were observed 
and no line splittings and shifts occur, a solution for "X. 
is obtained in the laboratory frame. 

In the calculations for mechanism (i) the mean- 
sq_uare values of the fiel-ds at the different nuclei, the 
correlation betv/een these fields and the veCLue of the 

correlation time^ are the unknovm factors in the absolute 
values of the elements of the relaxation matrix. For 
mechanism (ii) t is the only unknovm factor. 

G 

For mechanism (i), a fit with the observed spectra 
for strong irradiation can be obtained by performing the 
calculations for different values of f(i)/f(i'') (see Eq, 
III-4) and (see Bq. III-5). It kas been pointed out in 

the previous section tha,t the behavior of line 5 , in the 
double resonance spectra obtained under different conditions, 
indicates that the relaxation mechanism in the molecule does 
not connect symmetric and antisymmetric states. Ihe 
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relaxation meohanism? therefore^ seems to preserve the equi- 
valence of the two b' nuclei implied "by the spin hamiltoniaru 
Bor relaxation due to external random fields, the equivalence 
of the nuclei mil he maintained if the random fields at their 
sites are such that bS (t) = bS (t). This leads to the 

^1 "2 

condition that the random fields at and B^ are equal and 
completely correlated [I 6 ]. On the basis of this observation 
the calculations for mechanism (i) were performed by setting 
f(B^) = f(B 2 ) = f(B) and Cg 3 = Ij and varying only the two 

quantities f(A)/f(B) and C^^for high-field spectra"^. 

Bor mechanism (ii) the factors -|•(3GOs^*^JJ^>f-l) which 
determine the correlations between the different pairs of 
spins (see Eq. III-9) are calculated assuming the geometrical 
arrangement of the three spins in 2 , 6 -dibromoaniline to be 
the same as in an unsubsti tuted benzene molecule. The angle 
in benzene is 30 °, 30 ° and 120 °, giving 0 . 625 , 0,625 
and -0.125 for correlation between (ABg, B^B^), (AB 2 , 
and (AB^, AB 2 ) respectively. The distance between the protons 
is obtained (from the geometry of benzene) as 2,48 jf between 


+ It has been noticed in the computation that if f(B-|^) is 
set equal to f(B 2 ), a variation in ^5^32 P^odmee 

appreciable deviation in the intensity of line 3 from its 
single resonance intensity. The relative intensity and the 
line7/idths of other transitions in the double resonance 
spectrum are also insensitive to this constant as will be 
shown by a la.ter table (see Table V-5)o Thus even though 
*^BiB 2 equal to unity to conform to the theoretical 

requirements on the complete equivalence of the ‘'B^ nuclei, 
the experimental results do not allow an accurate estimate 
of this constant. 
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1 and B protons and 4.29 A* betY/een B protons. 

For every set of parameters describing tiie relaxa- 
tion process the matrix is obtained and substituted in 

Eq. (11-59) to obtain the peaJc signal intensity. All the 
other quantities in the numerator of that expression are 
obtained by diagonalizing ^ and calculating the matrix 

elements of I_^(i) and in that basis. It should be noted 
that the denominator in this expression contains the total 
linev/idth including contributions from inhomogeneity and 
instability as discussed earlier. This amounts to substi- 
tuting for in Eq. (11-59) just the value of the 

observed line width to obtain the value of • 

VD DISCUS SOT 

The calGUla.ted relative intensities in the double 
resonance spectrum for mechanism (i) with three different 
values 1.5,2^Q sued 2-5 for f (A)/f (B) and -C^=0.0 for all the 
three cases, and those for mechanism (ii) are given in 
Table (V- 2) along mth the observed relative inensities for 
the cases of irradia,tion of line 6 and 7 mth Vg = 2.2 Hz 
(see Fig- V-5) and for irradia.tion of line 1 with Vg = 5.25 
Hz (see Fig. V-7(a))- The value of given above is the 
final value chosen. Postponing the discussion on for 
the present, a comparison of the theoretical intensities 
for mechanism (i) with the observed values shows that 
f(a)/f(B) =2.0 gives the best fit with the experiment. 
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Tcille V-2 

Oise Lvecl rml oalculatod relative intensities for irradiation of 
linos 6 and 7 ?d.th. an amplitude of 2.2 Hz and of line 1 with an 
amplitude of 5*25 Hz. The random field calculations are for 
03^32 ~ = 0.0, For irradiation of linos 6 and 7, 

the intensities one expressed relative to the intensity of 
line 3 (S.F. 5) ac unity. 


Irradia- S.ht Observed 
tion at relative 

intensity 

tt 


Calculated relative intensity ± 


Random field 


f(l)/f(B) 


T 


1.5 

T 


2.0 


2.5 

“T“ 


Mpole- 

dipole 


5 


line 6 


V2=2.2 Hz 


1 

2 

5 

4 

5 

6 

7 

8 
9 

10 


0.33 

0.31 

0.30 

0.29 

0.28 

0.63 

0.67 

0.54 

0.62 

0.65 

0.42 

0.44 

0.42 

0.40 

0,40 

0.54 

0.53 

0.53 

0.53 

0.53 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

0.04 

0.04 

0.04 

0.06 

0.46 

0.44 

0.44 

0.44 

0.45 

0.20 

0.17 

0.20 

0.22 

0.19 

0.28 

0.26 

0.30 

0.33 

0.29 

0.09 

0.06 

0.07 

0.08 

0.07 

0.23 

0.11 

0.12 

0.14 

0.11 

0.37 

0.29 

0.34 

0.37 

0.36 

0.22 

0.17 

0.20 

0.22 

0.15 

0.55 

0.47 

0.46 

0.46 

0.51 

1.00 

1,00 

1.00 

1.00 

1.00 


0.04 

0.04 

0.04 

0.04 


line 7 
V2=2.2 Hz 


1 

2 

3 

4 

5 


6 
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Talle V-2(contd-) 


1 2 ' 

3 

4 

5 ' 

5 

7 

7 

0.55 

0.49 

0.50 

0.50 

0.51 

8 

0.33 

0.41 

0.39 

0.37 

0.41 

9 

0.50 

0.47 

0.45 

0.43 

0,50 

10 

0.20 

0.14 

0.14 

0.13 

0.14 

1 

- 

0.02 

0.02 

0.02 

0.02 

2 

0.26 

0.26 

0.27 

0.27 

0.28 

3 

1.03 

1.03 

1.03 

1.03 

0.92 

line 1 4 

0.88 

0,95 

0.97 

0.98 

1.20 

Vp = 5 . 25 Hz 5 

0 . 7 S 

0.66 

0.63 

0.62 

0.65 

6 

1.00 

1.00 

1.00 

H 

« 

O 

O 

1.00 

7 

0.25 

0.16 

0.16 

0.16 

0.23 

8 

— 

0.01 

0.01 

0.01 

0.00 


+ Serial numlDer of the double resonahce lines in a decreasing 
order of frequency (left to right in the spectrum). 

•|-^!rhese values are averages over several spectra (5--8) and are 
accurate approximately vfithin 5 per cent of the values given. 

^ Calculated relative intensities have been obtained by dividing 
the calculated integrated intensities by observed linewidths. 
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The calculated intensities for f(A)/f(B) = 1.5 and 2.5 do 
not shov; striking deviations from the observed values. 

However, it may he noted, for example^, for irradiation of 
line 6 with V 2 = 2,2 Hz (Table V-2), that the calculated 
intensities for f(A)/f(B) = 1.5 are consi stently higher than 
the observed ones for lines 1,2 and 3, and lower than the 
observed ones for lines 8,9 and 10, while the opposite is 
true for f(A)/f(B) = 2.5- Such a feature is noticed for 
many different cases of irradiation. This leads to the 
conclusion that f(A)/f(B) = 2.0 gives the best fit. The 
Values for mechanism (ii), however, also ^ree q_uite well 
with the experiment. Similar agreement is obtained for 
irradiation of these and other lines in the spectrum for 
various strengths of irradiation, so long as the strength 
of irradiation is large enough to satisfy Bloch approximation, 
indicating that the relative importance of these tv;o 
mechanisms for this molecule cannot be distinguished on the 
basis of this calculation, which depends only on the relative 
values of the elements of the relaxation matrix. 

The theoretical spectra for low values of Vg ^ 
however, depend on the absolute values of the relaxation 
terms. The observed and calculated relative intensities 
for irradiation of lines 6 and 7 with Vg = 0,1 Hz are given 
in Table (V-3). The calculated relative intensities for 
mechanism (i) are for f(A)/f(B) = 2,0, = 0.0 and three 

different values 0.1, 0.2 and 0,3 rad./sec. for TfCB)!: 
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values are averages of alout five spectra said are accurate approximately within 5 per cent 
of the values given. ^ r ^ 






128 


(this was the oommon factor not included in the calculation 

of relxation matrix elements for spectra vdth large Vg) and 

those for mechanism (ii) axe for t = 5.0 x 10“^*^ sec. and 
-11 + 

2.7 X 10 sec.^ Prom the trend of the relative intensities 

given for mechanism (i) and hy mahrirg a similar comparison 

for irradiation of other lines, it was found that the spectra 

calculated with 2f(B)T: =0.2 ra„d./sec. show consistently 

c 

good agreement. 

The results for mechanism (ii) with t = 5.0 x 10"^*^ 

o 

sec. also agree well with the experiment. The correlation 
time appropriate for this mechanism is that for rotational 
diffusion , for which a hard sphere model for the rotation of 
the molecule gives! 


^ _ 4-TT--n 

^c - 3hT 


(V-1) 


t These calculations \7ere performed by assuming tha.t only 
one off-diagonal element of % between the states 
connected by irradiation, is non-zero (see Sec. II0-2(a))- 
However for irradiation of lines in ^ B ^ group in which each 
line has a close neighboiir, more than one off-diagonal 
element of yC could become significant. To estimate the 
effect of this a calculation in the rotating -frame 
assuming two off-diagonal elements of % , between the 
states giving the close doublet (the doublet which is being 
irradiated), to be Bimultoneously cignificant was per- 
formed (see Sec. IIC-2(b)). This calculation leads to 
values of X which axe not significantly different from 
those with only one off-diagonal element. 
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where rj is the viscositj oJid '' the hard sphere radius [?]. 

For the present sojnple p -ns 1,0 oentipoise. At room tempe- 

-1 0 

rature, the value of t = 5,0 x 10 sec, leads to a. hard 
sphere radius of about 8,0 A*. This value of is consi- 
derably larger than the molecular dimensions. Furthermore, 
the natural linewidth of line 3? for example, given by 
R . //tc (full width at half height in Hz) becomes 0.65 Hz, 
for the above value of x , which is larger than any of the 

G 

observed linevhLdths, including the inhomogeneity, for single 
resonance transitions. On tho other hand, if one assumes a 
value of ''a’' of the order of the molecular dimensions 
(2 - 4 a''), the calculated intensities (see fable V-3, for 

T = 2,7 X 10“^^ sec. corresponding to ^a"^ eq^ual to 3,0 A) 
do not agree mth the experiment, These considerations show 
that mechanism (ii) is not a dominant reladation process for 
this molecule. It should^ however^ be noted that the above 
argument is based on the estimn.tes of correlation time for 
the hard sphere model of rotational diffusion, and therefore 
the Values given should be considered yd. thin the limitations 
of this model. 


The effect of varying the correlation constant 0^ 
on the calculated spectrum for mechanism (i) is given in 
fable (V-4) along y/ith the natural linewidthS(R^^^^^,/a) of 
the different transitions for different values of Also 

shown, in the fable (¥-4) are the observed relative intensities 
and the sum of relaxation and instability widths , obtained by 
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lable V-4 

Otserved and oalcnlated relative intensities (pealr-heights) and 
line?7idtiis of tiie doifole resonance transitions when line 6 is 
irradiated with an amplitude of 2»2 Tiz. I'he caloulated inten- 
sities and linewidths are for mechanism (i) having 0 -d -□ = 1.0, 
f(A)/f(B)=2.0 and 2f(B)T^=0.2 rad./sec. ^1^2 


S.l 


Intensities 



linewidthst 


it Ohser 
ved 

Caloulated 

Oh served"^ 

Calculated 
relax, width 





^AB 


1.0 

0.5 

0.0 

Total 

Eel ax . + 
instah . 

1.0 

0.5 

0.0 

1 

0.33 

CM 

• 

O 

0.30 

0.30 

0,65 

0.29 

0.35 

0.31 

0.27 

2 

0,63 

0.61 

0.63 

0.64 

0.78 

0.42 

0.35 

0.30 

0.25 

y 

' 0.42 

0.40 

0.41 

0.42 

0.83 

0.47 

0.34 

0.29 

0.24 

4 

0.54 

0.52 

0.52 

0.53 

0.54 

0.36 

0.24 

0.25 

0.25 

5 

1.00 

1.00 

1.00 

1.00 

0.60 

0.42 

0.13 

0.13 

0.13 

6 

- 

0.04 

0.04 

0.04 


- 

0.24 

0.25 

0.25 

7 

0.46 

0.44 

0.44 

0.44 

0.62 

0.44 

0.26 

0.26 

0.25 

8 

0.20 

0.20 

0.20 

0.20 

0.50 

0.50 

0.15 

0.19 

0.24 

9 

0.28 

0.29 

0.29 

0.30 

0.42 

CM 

O 

0.17 

0.21 

0.25 

10 

0.09 

0.07 

0.07 

0.07 

- 

- 

0.19 

0.23 

0.27 

t 

Serial nounhers, 

as in 

Tahle 

(T-2). 






t Bull linewidths at half height in Hz. 

fi- These values are averages over several spectra (5-8) and 
are accurate approximately within 5 per cent. 
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sutstracting tlie iniiomogeneity contribution from the observed 
widths. The intensities are not sufficiently sensitive to 
the velue of but the values of this 

constant quite considerably. It may be seen from the 
calculated values of ~ 

that the effect of grea,ter correlation between A and B is to 
introduce an as 3 rnimetry in the natural linewidths in the 
spectrum. The observed relaxation plus instability' Widths, 

however, do not show any systematic variations and since the 

■ffsa 

instability is assumed to contribute same vadth to all the 

A. 

transitions a comparison of the observed and calculated widths 
shows that is close to zero. Tliis estimate is only 

approximate as it is coupled vnth the estimates of inhomo- 
geneity contribution made from those linewidths. The effect 
of variation of 0-^ -n is shovm in Table (Y-5). This table 

V2 

shows that neither the calculated double resonance relative 
intensities nor the calculated linev/idths are sufficiently 
sensitive to this parameter. However, as remarked earlier 
C-n -n has been token as unity in the present case to maintain 

V2 

the complete equivalence of the two B protons. ^ It may be 
noted that the values of correla,ticn constants = 0 and 
Otj -n =1.0 have been finally deduced in this work by 

V2 

experimental evidence which is not striking enough to make 
a sharp dis tine tiai between the different values* The values 
given may , therefore , be considered to be indicative of the 
degree of correlation, rather than to signify an accurate 


estimate of the same. 
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lable V-5 

Calculated relative intensities (peal:™ li eights) and 
relazation v;idths of the doulle resonanoe transitions 
when line 6 is irradiated hy an amplitude of 2.2 Hz. 
Ihese calculations are for mechanism (i) having 
f(A)/f(B) = 2.0, = 0.0, 2f(B)T = 0.2 rad./sec. 

and different valu^ of C-^ , 

BlBg 


S.l. 

Intensities 

linewi dths 





“BiBj 


1.0 

0.5 

0.0 

1.0 

0.5 

0.0 

1 

0.30 

0.30 

0.31 

0.27 

0.24 

0.21 

2 

0.64 

0.65 

0.65 

0.25 

0.25 

0.24 

3 

0.42 

0.42 

0.42 

0.24 

0.22 

0.20 


0.53 

0.53 

0.53 

0.25 

0.22 

0.19 

5 

1.00 

1.00 

1.00 

H 

• 

O 

0.16 

0.19 

6 

0.04 

0.04 

0.04 

0.25 

0.24 

0.22 

7 

0.44 

0.44 

0.44 

0.25 

0.22 

0.19 

8 

0.20 

0.20 

0.20 

0.24 

0.22 

0.20 

9 

0.30 

0.29 

0.29 

0.25 

0.25 

0.24 

10 

0.07 

0.07 

0.07 

0.27 

0.24 

0.21 





133 


The -theoretical values of for the c?.se of 

nechaJiism (i) having f(A)/f(3) = 2,0, 2f(B)TQ = 0.2 rad. /sec. 

C.-D = 0.0 and C-n t. =1.0 are plotted in jiixtaposition with 

Ah 2 

the double resonance spectra in Figs. ('V-3) - (V--7). There 
is good agreement between the theoretical and experimental 
spectra in all the cases. 


It shoul'd be pointed out that if the effect of 
magnetic field inhomogeneity is not properly taken into 
a,GGOunt the values of the relaxa-tion parameters determined 
and the possible distinction betv/een the theoretical calcu- 
la,tions for different mechanisms would be completely altered 
and may lead to erroneous results. For example, in the 
present case^if the magnetic field inhomogeneity is assumed 
to contribute the some v/idth to the different double resonance 
transitions, the value f(A)/f(B) may have to bo changed to 
nearly 4.0 and to 0.5 to fit the obseiwed spectra. 

It must be emphasized that it has been possible to 
conclude that mechanism (ii) is not an important relaxation 
process for the present case, only by combining the analyses 
of the results of double resonance spectra obtained for both 
strong and ?ireak irradiation. The spectra for strong irra- 
diation contain, in general, more information in the sense 
that there are more transitions for ?\rhiGh a fit con be 
obtained, but the final fit in this case depends only on the 
relative values of the elements of the relaxation matrix. 

The spectra for weak irradiation, on the other hand, contain 



134 - 


fewer parameters to obtain a fit, but the calculation depends 
on the absolute values of the elements of the relaxation 
matrix, This is an interesting distinction between these 
two categories of spectra and can be used vd.th considerable 
advantage in the study of relaxation by double resonance 
methods. 

This analysis loads to the conclusion that the 
random isotropic fields at the site of nucleus A have a mean 
square value twice thai at the sites of B nuclei. This 
excludes the other spins in the molecule as important sources 
of these random fields as nucleus A is fa.rther off than the 
B nuclei from the remaining spins in the molecule. The 
proton relaxation thus appears to be caused primarily by 
intermole cular dipolar interactions. Since the sample is a 
20 per cent solution by weight in CCl^ (wliich is nearly a 
saturated solution) en.ch molecule of 2, 6-dibromoaniline is, 
on- an average, surrounded much more by CGl^ molecules than 
the molecules of the same kind (a 20 per cent solution by 
v;eight of 2, 6--dibromooniline in CGl^ has approximately seven 
GGl^ molecules per one 2, G-dibromo aniline molecule). It 
was not feasible to look for a, possible conce'^ntration 
dependence of this relaxation since the range in which the 
concentration could be reduced, mthout seriously affecting 
the signal to noise ratio, would slightly increase the 
number of GGl^ molecules and would not significantly alter 
the physical situation in the immediate vicinity of a 
dibromo aniline molecule. 
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VI - DOUBLE EESOIillCS SIUIY 0 ? PROIOH HBIAKAIIQH IH A I¥0 
SPM SYSTEM (AB) COUPLED 10 A QUADRIXPOLAR lUGLEUS 

VIA nilRODUGIION 

In this chapter the restilts of a proton-proton douhle 
resonance study of relaxation in a two spin system (AB), the 
spectra of which exhibit broadening due to coupling with a 
neighbouring nucleus, are presented. The AB system is 
formed by the protons in 2-bromothi azole which is studied in 
the form of a neat liquid. The scalar coupling of the protons 
with the nucleus modulated by the rapid quadrupolar rela- 
xation of the latter is an important source of proton relaxa- 
tion in this molecule. The analysis of double resonance 
experiments performed by irradiating each of the single 
resonance transitions for different strengths of irradiation 
led to some detailed information on the scalar coupling 
interaction viz. the magnitudes and relative sign of proton- 
U coupling constants, as well as the parameters describing 
the other mechanisms of relaxation in this molecule, in 
particular internal dipole-dipole interaction. Conventional 
relaxation time and linewidth measurements were also used to 
facilitate the double resonance analysis. It was shown by this 
study that, in general, it is possible through relaxation 
studies by double resonance, to obtain the relative sign of 
the coupling constants of protons with quadrupolar nucleus. 

The spectra for low strengths of irradiation were analysed 
by the graphical method, mentioned in Sec. IIIC. 
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VIB EELAIITE SIGH OF SCjaAR COUPIUTG CONSIMIS 


When there is more than one spin coupled to a q_uadru- 
polar nucleus, the relaxation matrix elements for the spin 
system contain terms arising from cross-products hetv/een 
various scsiar coupling interactions. These interactions are 
completely correlated because the time- dependence is caused 
by the same quadrupolar nucleus and these cross-terms do not 
therefore vanish. For two spins both of which are coupled to 
the same quadrupolar nucleus the scalar coupling interaction 
may be written as (see Eq_. III-IO), 


UK) + Jjj, I(B)>3,(t) 

( VI-1 ) 

The spectral densities (Eq_. 11-30) will then contain terms 
of the type; 






1(A) + 


I(B))-Sjj(t) 


» Sjj ( t-T ) 




{VI-2) 


which depend^in general, on the relative sign of and 
This suggests the possibility of determination of these rela- 
tive signs by studying the relaxation matrix elements. 


These relaxation matrix elements depend on the relative 
sign of and even in the case of single resonance. In 
single resonance spectra these matrix elements determine the 
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linevddtiis of the different transitions. Hov/ever , in cases 
in v/hich ^nd ( - OJ'g)^T^'^l the 

relazation is oontrihuted only hy the terms 2'n:[j^I^(A) + 
JBj;fl2(B)]s|(t) of Bo. (VI-1) (see Sec. IIIO) and the widths 
of each of the transitions belonging in the zeroth order, 
to V nucleus is 4"^^ (J^-rCOs^O + JBjjSin^©) and to the 'B'' nucleus 
is 471;^ (JbjjCos^© + J^jsin^e) where tan © = 

Ihe single resonance spectra in such cases do not, therefore , 
depend on the relative sign of and 

The double resonance spectra, however, depend on a 
variety of relaxation matrix elements and are therefore 
expected to depend on the relative signs through the intensi- 
ties as well as linemdths of various transitions. This fact 
has been utilized in the present case to find the relative 
sign O-j^ J and • 

VIC EXPERBffihTAL EESDITS 

1 . Single Resonanc e S£ e c trum 

The single resonance proton spectrum of 2-bromo- 
thiazole is shorn in Eig, (IV-5). This is a typical four 
line AB spectrum, having additional broadening due to scalar 
coupling of the protons with ni'crogen. The low field ('a" 
proton) doublet is assigned to the proton in position 4. 

This proton being nearer to nitrogen compared to the proton 
in position 5 (''B' proton) will couple strongly with it, 
giving a broader doublet. This assignment is same as that 

in [44] . 
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The spectral analysis of an -AB spin system is straight 
Lorward [2]. The analysis gives the values of the coupling 


constant = 3.55 + 0.05 Hz and the chemical shift 
= 18.8 + 0.1 Hz between the two protons. 


A B 


2. T^ and T^ Measurements 

Measurements of spin-lattice relazation time T^ and 
spin- spin relaxation time T^ were made in order to derive 
some information on the magnitudes of scalar coupling of the 
protons with . The magnetic field v/as made sufficiently 
inhomogeneous for the doublet of eanh proton to collapse into 
a broad line and T^ and Tg for each such doublets were 
measured. T^ was measured by the adiabatic rapid passage 
method [45]. The result of one such experiment is shown, in 

Big. (VI-1) and a plot between log (M -M ) and time in 

^ z 

Big, (VI-2). The slope of this straight line is a measure of 
T^. The Values of T^ obtained by tliis method for the * A^ 
proton (lov\r field broad doublet) and the ' B proton (high 
field sharp doublet) are 9.1 + 0.5 sec. and 8.8 + 0.5 sec. 
respectively.. 


The spin- spin relaxation time T^ wa^s measured for the 
^B*^ proton by adiabatic rapid passage into the center of a 
u-mode resonance with a large r.f. field (large compared to 
the inhomogeneity of the magnetic field), stopping the sweep 
at that point and monitoring the rate of decay [45]. This 
decay is exponential with a tine constant T^ and a typical 
curve obtained for the B doublet is shovm in Big. (VI-3), 
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VI-1. 


Re o every oiirve lor liie adi.r'batic 
F.e1;hod of neroureiii nt for the 
douhlet 01 2-1 rorac till c;3 ole. Ilie 
peahs are one oco ond aparL', 


rapid passage 
high field 
suoocssive 
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Pig. VI-2. Iog(M - M ) -versus time (t) graph for I-, 

Q 3 J- 

meo-suremcp.t of high field doublet of 
2-b r omo thi az ol e . 
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Fig. VI- 3. Saturation curve for the adiabatic rapid 
passage method of 1 2 measurement for the 
high field doublet of 2-hromothi azole . 




t — 

IN SEC 
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from ’■.viiicii a value of is obtained as 3.6 + 0.2 sec. The 

value of Tg proton is conpar actively very small 

14 

due to its stronger coupling v/ith IT and therefore could 
not he measured hy this experiment. 

nitrogen resonance was observed on a Varian 

V-4200 wi deline n.n.r. spectrometer, operating at 2.695 

(field -^8, 16 K.G.), using audio sweep modxilation. Modulation 

frequency of 20 Hz and sweep amplitude of 1.0 gauss were used. 

She spectrum was observed in the v-mode and the peah-to-peah 

linevrf-dth of the first derivative, after applying a modulation 

correction [46], was obtained as 0.98 gauss, which gives a 

fullwidth at half height of a lorentzian as 1.70 gauss. 

This leads to the relaxation time t (T-, = Tp = t ) for 

q J. c. q 

as 0.61 millisec. 

3. Double Resonance 

Some of the typical double resonance spectra obtained 
are shown in Digs. (■V'I-4)-( VI--6) . The calibration for the 
strengths of irradiation was done in the manner described in 
Sec. YC. Figure (YI-4) shows the irradiation of each one 
of the single resonance lines by an amplitude of 0.135 Hz. 
'low"^ irradiation amplitude spectra have been obtained for 
two more strengths viz. 0.076 and 0,043 Hz, for irradiation 
of each of these lines. These spectra show significant 
intensity changes in the spectrum of the proton other than 
the one being irradiated. These changes are large for transi- 
tions having energy levels in common with the pair being 
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Eig, ■VI-4. Frequency sweep double resonance spectra of 
2“bromot]aia2ole, for '‘'lov;^* strengths of 
irradiation, obtained by irradiating (a) inner 
line of proton a, (line '5)s (1) onter line of 
proton A, (lino 4), (c) inner line of proton 
B,(line 2) and (d) outer line of proton B, 
(line 1). file strength of irradiation in 
each case is 0,135 Hz. 



8 
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Fig. VI-5. Observed ajid ooloiilated double resonance 

spectra for irradiation of (a) outer line of 
proton A (line 4) and (b) inner line of proton 
A (line 3) in 2-bromotbi azole, v»dth amplitude 
of iriadiation 2.4 Hz. Tbo calculated spectra 
are for internal dipole-dipole interaction 
and scalar coupling (mechanism (i)) with 
Iv — 4.0 and jj'i ~ 0.3. 




o 
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Fig. VI-5. Observed and oalcnla-tod double rosoncnco 

spectra of 2-bromotki azoic for irradiation 
of (a) outer line of proton A (line 4) vri.th 
strength of irradiation 1.6 Ha and (b) outer 
line of proton B (line 1) mtii strength of 
irradiation . 1,75 Ha. The calculated spectra 
are for internal dipole-dipole interaction 
and scalar coupling (mechcnism (i)) with 
K = 4-.0 and = 0.5. 
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irradiated and can be qualitatively understood on the basis 
of argunentSj, using an energy level diagram, similar to the 
ones described earlier (Sec. VC-1), She spectra v,d.th ' high^ 

strengths of irradiation shov/ all the features of a double 

* 

resonance spectrum (Figs. VI-5,6), She contribution to the 
linevd-dths to these double resonance transitions is also, in 
general, different for different transitions. Eeliable 
estimates of this contribution for different transitions are 
obtained by using the plots of versus . These 

estimates are, hov/ever, approximate since the spin system , 
in the present case, is not vjealcly coupled. 

VIE miJjYSIS 

1 • Eatimates of Scalar Gounling Strengths 

The contribution of either dipole-dipole interaction 
or external random fields to T^ and T^, for very short corre- 
lation times, are equal (extreme narrowing), Virhile this is 
usually not the case for the scalar coupling mechanism, since 
the corresponding correlation time is rather long in 
comparison [?]. Therefore^ wrherever scalar coupling contri- 
butes to relaxation, a difference between the values of T^ 
and T 2 gives information on the extent of this coupling. 
Estimates of coupling constants and of protons A and 

B respectively with may then be obtained by using Eqe. 
(111-49,50) and an assumption that the relaxation matrix 
element of the type represents where a p 
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is a '' transition in the zeroth order"^, ihis leads, for 

a coupled spin systen AB under the assumption that only terms 

of the tjTpe I„S of .(t) in Eq_. (VI-1) are significant, 

Z Z t s c 

to: 


(1 _ 1 ) 
^2 


and 


[ik )o.3 


^2 B 


+ J|j, oos=0) 
(VI-5) 




,1 


An^Si S+1) 


(J^ - J^y) cos 20 ^ (VI-4) 


where the subscripts ''a” and *“ B refer to the two protons and 
the subscript "sc'' to scalar coupling. Equation (VI-4) 
gi\’‘es the difference between the observed single resonance 
linewidths of A and B protons, if the contribution to these 
linev/idths from other sources (namely, inhomogeneity of the 
magnetic field and other mechanisms of relaxation) is assumed 
to be the same, 

f It may be noted that it is, in principle, possible to 

derive the values of the scalar coupling constants by double 
resonance data alone. However, estimates made on the basis 
of approximate determinations of and I 2 serve as a 
useful starting point for obtaining fit with the double 
resonance experiments, and v/ould facilitate the determina- 
tion of the parameters of the other mechanisms. 
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Using the experimentally deternined values of 
and (‘l 2 ^B respectively 8.8 sec. azid 3*6 sec. and = 0,61 
nillisec, in Eq. (VI-3), end equating the observed linewidth 
difference hetween the A and 3 proton resonances 0.35 Hz) 
to Eq. (VI-4), gives: = 8.8 Hz and = 3*1 Hz. Ihis 

value of compares favourably \vith a similc.r estimate of 

10.8 Hz by Kintzinger and lehn for thiazole [37]. 

For these values of and and the observed t^, 

the condition mentioned in Sec. IIIC, for complete washing 
out of the scalar coupling 1) is satisfied, ThereforSjfor 

the present sample the treatment o_f scalar coupling only as 
a. mechanism for relaxation, is valid. 

2. Other Details 

Ihe proton-proton distance in this molecule enters 
the internal dipole-dipole intera.ction calculation. Informa.- 
tion is available on bond angles and bond lengths in thiazole, 
through studies of C^^-H coupling constants [47]. Assuming 
that the bromine substitution does not significantly alter 
the geometry of the molecule, the distance between the two 
protons in 2-bromothi azole is obtaine-d a-s 1.66 . 

Theoretical single and double resonance spectra have 
been calcilated by considering (i) internal dipole-dipole 
interaction and scalar coupling and (ii) external isotropic 
random fields and scalar coupling as possible alternative 
mechanisms for relaxation. For double resonance spectra with 
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’'high'' strengths of irradiation, for which jYH2{ j ["aa' |’ 

Bloch approximation has been used (see Sec. IIC), Since 
under Bloch a,pproximation the calculated spectra depend only 
on the relative value of the relaxation matrix elements, 
these spectra have been calculated by varying a parameter K, 
v/hich determines the relative contribution of dipole-dipole 
and scalar coupling in mechanism (i) and random field and 
scalar coupling in mechanism (ii). Bor random fields, however, 
there are two more parameter namely the ratio of moan 
square fields f(A)/f(B) at the sites of the two protons and 
0 the correlation constant between these fields. The ratio 
f(A)/f(B) was varied in the c alciilations to obtain a fit with 
experiment, but C was fixed from the information available 
from the observed single resonance linewidths. The observed 
single resonance linewidths are more for the tvro inner lines 
in the AB quartet, in this case, than the outer lines, indi- 
ca,ting that C is close to unity [ref. 7, p. 509], and was 
takien to be 1. It may further be remarked in this connection 
that while the internal dipole- dipole also contributes greater 
vd-dths to the inner lines compared to outer ones, in the AB 
quartet, the scalar coupling contributes equal widths to the 
two lines of a proton resonance in it. 

« // 

Bor the double resonance spectra vhth low strengths 
of irradiation, the observed and calculated straight line 
plots (Eq, 11-80) were drawn and compared. 
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VIE Discussion 

■Table (VI-1) shows the observed and calculated double 
resonance peah heights, the observed single and double reso- 
nance linevri-dths and calculated single and double resonance 
relaxation mdths for irradiation of outer line of reso- 

nance by anplitudes of irradiation Vg = 1.6 and 2.4 Hz and 
for irradia-tion of outer line of B resonance by an ampli- 
tude V 2 = 1.75 Hz. A satisf antory fit for observed peal 
heights and linewidths is obtained for nechanisn (i) for K = 
4.0 and = 0.5 and for nechanisn (ii) for K = 4.0, 

f(A)/f(B) =0,5. The parameter K is equal 

to for mechanism (i) and equal 

^ i ■ 

to [ (le-re^ g/3)/( 2f (B)v^) for mechanism (ii). A variation 

in does not seriously affect the calculated peak- 

heights, but alters the calculated linewidths. The 

versus A. plots for approximate estimation of inhomogeneity 
J 

contribution to various double resonance transitions were 
dravm and one such plot for irradiation of ''A'' proton by Vg = 
2,4 Hz is shOYCi in Pig, (71-7). This plot indie ates, for 
irradiation of the outer line of '’A'" resonance which corres- 
ponds to A^ =-1,94 Hz, approximately 40 per cent increase in 
the recoverable part of single resonance inhomogeneity 
contribution (the instability contribution is not recoverable) 

Though this quantity has been estimated from T^ and T^ 
measurements, the accuracy of these measurements is limited 
in as much as these measurements were made for overlapping 
doublets. 
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ITig. VI~7. Plots of caloiilated double resonance frequen- 
cies of ' B transitions versus " frequency- 
offset of irradiation of nuclei. These 
plots are calculated for the AB spin system 
for Vg = 2.4 Hz. Irradiation of line 1,2,3 
and 4 correspond respectively to = 20.74? 


17.19 ,1.61 and -1.94 Ez 
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for- lines 1 and 2 of "the double resonance spectrum and a 
corresponding decrea.se in lines 3 and 4. If the recoverable 
pant of the inhomogeneity is assumed to contribute a full- 
vri-dth at half height 0.15 Hz to single resonance transi- 
tions (see Sec. 70-2), the observed and calculated single 
and double resonance vddths fit best for = 0,3* A 

detailed comparison of the double resonance spectra obtained 
for different frequencies and amplitudes of irradiation 
beon - m - a -de and led to a satisfactory fit for the above mentioned 
parameters for mechanisms (i) end (ii). 


Ihe effect of a change of sign of with respect 
to is also shown in fable (7T-1). The single resonance 
linewidths are not changed as shown in Sec, (VIB). However, 
the double resonance linewidths sho?; a marked difference 
between the two alternative cases of s.ame sign and opposite 
sign. The linemdths for the case of opposite sign are widely 
different from the observed ones, leading to a conclusion that 
both and Jgjjj- are of the same sign. 


The spectra for *low amplitudes of irradiation have 

been calculated for the above two mechanisms with parameters 

ihose. 

K and same as^which led to a satisfactory fit for 

high" amplitude of irradiation case, and by varying the 
absolute values of the relaxation matrix elements. This 
amounts to, for a given value of a simultaneous variation 
in and for mechanism (i) and in 2f(B)'i;^, and 

mechanism (ii) such tha,t K and u^^e unchanged. 

A fit for '‘low' and 'high' amplitude of irradiation spectra 
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then determinos tiie "vaiues of these parameters as ottained 

by double resonance analysis. Hcv/ever, sonc infoniia.tion on 

and is available from rel.axa,ticn tine measurements 

and therefore the absolute values of the relazaticn matrix 

elements need be varied only over a small range. Table 

(VI--2) gives the observed and calculatod values of R, 

Eq. (11-79), for various transitions in the double resonance 

speotra for ' lov; ^ irradiation of A group of lines. The 

calculated values ane for Z = 4.0; = 0.3 and 

different absolute values of relaxation matrix elements. A 

satisfactory fit vdth the observed values of R is obtained 

for mechanism (i) mth t = 7.0 x lO""^^ sec. Ror mechanism 

o 

(ii) no satisfactory fit is obtained. Also shov/n in Rig, 
(VI-8) are the plots of observed and calculated ratios 
(R+1)/(R-1) between various transitions versus l/(vp, for 
irradiation of the ' group of lines vdth ''low^ amplitudes 
of irradiation (see Eq_. 11-80). While the observed ratios 
are indicated by various points on the graph, the calculated 
ratios are shown by straight lines. The fit betvreen the 
observed points and the calculated radios for mechanism (i) 
(continuous lines), for K = 4.0, = 0.3 and t = 

7.0 X 10 sec. is satisfactory. These parameters give 
^ '^q_ millisec. Ror 

mechanism (ii), however, no satisfactory fit could be 
obtained for any plausible values of the parameters, and the 
plots in Rig. (VI-8) for this mechanism (broken lines), 
for K = 4.0; = 0.3 and 2f(B)T;^ = 0.19 rad./sec. 
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lalDle VI- 2 

^observed and C8.1cv.lated values of E (see Eq .. 11 - 79 ) for 
'low" irradiation of 'A' group of lines, IIig calculated values 
for- Tooth. rGchanisms (i) and ( ii ) are for K = 4.0 and IbI ' T / JaU = 
0 . 3 - Eor mechanism (ii) the additional parameters are f(A7/f(B) 
= 0.5 and C = 1 . 0 . 


Irrad- 

iation 

3 , 1 !; 



Ratio 

of 

linos 


Oh so 
rvod 


mechanism (i) 

___ 

X 10 SGC 

c 

11.0 7.0 4.7 


mechanism (ii) 
2f(B)v^ 

0.29 0.19 0.15 


2/1 

0.043 

^ 

2/1 

line 4 0,076 

3/1 

2/1 

0.135 

3/1 

2/1 

0.043 

^ 

2/1 

Line 3 0.076 

4/1 

2/1 

0.135 

4/1 


0.76 0.88 0.77 0.64 0.94 0.'89 0.81 

0.89 0.94 0.88 0.82 0.93 0.87 0.77 

0,61 0.72 0.59 0.48 0.86 0.76 0.68 

0.79 0.86 0.79 0.73 0.83 0.73 0,63 


50 

0.54 

0.45 

0.40 

0.73 

0.64 

0.-59 

71 

0.76 

0.71 

0.69 

0.69 

0.59 

0.53 

33 

las 

1.36 

1.67 

1.06 

1.12 

1.23 

24 

1.10 

1.20 

1.38 

0.98 

0.96 

0.92 

89 

1.47 

1.84 

2.30 

1.16 

1.28 

1.43 

48 

1.27 

1.48 

1.74 

0.94 

0.90 

0.84 


2,54 2.04 2.50 2.88 1.33 1.48 1.60 

1.89 1.59 1.86 2.07 0,88 0.81 0.78 
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Fig, ¥1-8. 


Plots of observed oJid c. 3 loifL oted ratios (R+1)/ 

(E-l) versus l/v^ s^3?engths of irradia,- 

tion, for 2-bro3iothiaaolG . The poincs indicate 

the observed ratios and the straight lines the 

theoretical ratios. The solid lines are for 

mechanism (i) mth K = 4.0, = 0,3 pn.d 

-1 ? 

X = 7.0 X 10 sec. and the broken lines are 
c 

for mechaJiism (ii) vd-th K = 4.0, 0-3j 

2f(B)i: = 0.19 rad. /sec., f(A)/f(B) = 0.5 and 

G 

0 ,-n = 1.0. For irradiation of line 3 A A 

Ad 


represents, in the definition of S, the ratio 
of line 4 to line 1 and #8 the ra,tio of line 
2 to line 1. Tor irradiation of line 4 OC 
represents the ratio of lino 2 to line 1 and 
A B the ratio of line 3 to lino 1. Foi^ each 
place of irradiation three strengths 0.043, 
0.076 and 0.135 Hz were used. 



4- 1 } /( K 
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( ^or- 9.7 Hz ojid 2,9 Hz) show large deTriations 

from experinertal values. These considerations show that 
the external random fields along wath scalar coupling do not 
account for the chserved cliaagos in the spectra with low'"" 
amplitudes of irradiation and hence nechanisn (ii) is ruled 
out. The internal dipole-dipole interaction along vdth 
the scalar coupling accounts for all the features of the 
single and double resonance spectra, and is , therefore, the 
dominant nechanisn for relaxation in this molecule. As the 
relaxation contribution of internal dipole-dipole interaction 
is strongly dependent on the distance betv/een the interac- 
ting spins (r“^), this result is not surprising in view of 
the snail proton-proton distance (1,66 A®), in this molecule, 
(compa,red to 2.48 betv/een the ortho protons in benzene). 
The viscosity of 2-bromothi azole was measured at room tempe- 
rature and a value ^2,0 centipoise obtained. TMs leads 
to a hard sphere radius (see En. Y-1) of the order of 2 II , 
which fits well vhth the molecular dimensions in this case. 

The theoretical spectra for mechanism (i) having 
K = 4.0 and = 0.5 are plotted in jinctaposition with 

the double resonance spectra in Figs. (YI-5 and 6). These 
show good agreement betv/een the observed and calculated 
spectra. 

Table (YI-3) gives the contribution of internal 
dipole-dipole interaction and scalar coupling to the single 
resonance relaxation matrix elements of the type 


and 
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Table VI- 3 

Contributions of internal dipole-dipole interaction 
and scalar coupling to various single resonance 
relajxation matrix elements (in rad. /sec. ) for mecha- 
nisra (i) having K = 4.0, JgTvj-Aj^j = 0.3 and = 

7.0 X 10“^^ sec. 


Matrix 

elcnentt 

Dipole- 

dipole 

Scalar 

coupling 

Total 

12 12 

-0.174 

-0.149 

-0.323 

13 13 

-0.149 

-1 . 508 

-1.657 

14 14 

-0.171 

-2.565 

-2.737 

2 3 2 3 

-0.076 

-0.731 

-0.808 

2 4 2 4 

-0.174 

-1.508 

-1,681 

3 4 3 4 

-0.149 

-0.149 

-0.298 

1111 

0.171 

0 . 000 

0.171 

112 2 

0.035 

0.000 

0.035 

113 3 

0.021 

G.OOO 

0.021 

114 4 

0.114 

0.000 

0.114 

2 2 2 2 

-0.086 

-0.013 

-0.099 

2 2 3 3 

0.018 

0.013 

0.031 

2 2 4 4 

0.034 

0.000 

0.034 

3 3 3 3 

-0.065 

-0.013 

-0.078 

3 3 4 4 

0.023 

0.000 

0,023 

4 4 4 4 

-0.171 

0.000 

-0.171 


The numbering of the states is such that for the 
i\3 spin system, in the zeroth order these states 
are given, in the increasing ordej; by, aa, ap, Pa, 
and pp. 
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parsneters for ^/rhich a soJtisf actcry fit of 
the double resonance spectra has bean obtained. This table 
shows the relative contribution of these two mechanisms to 
the relaxation in this spin system. 

It is concluded from this study that the internal 
dipole-dipole interaction and the scalar coupling of the 
protons mth nitrogen, which is relaxing rapidly due to 
quadrupole relaxation, in 2-bromothi azole, accomtsfor proton 
relaxation in this molecule and that the external randon 
fields do not contribute significantly to the relaxation 
processes. Various parameters describing the internal 
dipole-dipole interaction (t 5=3 7.0 x 10“^^ sec.), and the 
scalar coupling mechanism 9.7 Hz, p 2.9 Hz and 

of the same sign) have also been obtained. 
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SUiaiARY iND GOIIOLUSIOL'S 


Relaxation processes in coupled opin systems have 
"been- studied by steady-state homonuclear double resonance 
experiments. The systematics of the double resonance method 
and details of the relaxation processes have been examined. 
Double resonance experiments have been performed on a coupled 
symmetrical three spin system (ABp) and on an AB spin system. 

It is foimd that it is useful to study double 
resonance spectra obtained over a wide range of strengths of 
irradiation, as the solution of the density-matrix equations 
obtained for conditions of high. r?nd lov/ irradiation 
amplitudes possess some distinctive features which may then 
be used with considerable advantage to distinguish bet?/een 
different relaxation mechanisms. A detailed consideration 
of the inhomogeneity of the magnetic field is sh07«i to be 
of crucial importance for the analysis of the relaxation 
effects and a m-ethod of including it in the density matrix 
calculations of the double resonance spectra is suggested, 

’The inhomogeneity of the magnetic field is found to contri- 
bute a full -width at half-height of the order of 0,15 Hz to 
single resonance transitions^ in the present experiments, which 
was almost wholly recoverable in double resonance spectra 
under suitable experimental conditions. 

Isotropic random field model, convenient for descrip- 
tion of relaxation effects in such studies, is shown to be 
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a plausilsle model for -various intermolecular and intramole- 
cular interaotioiis under conditions normally realised in 
experiments. Wlienever such e, model is applicable the cons- 
tants describing the correlation between the random fields 
at the sites of different spins have been discussed. A 
definition of this constant jV/lii oh includes the possibility 
of completely correlated fields ■v'n.th different mean square 
values, is used. 

Ihe relaxation of the ABg spin system of tho ring 
protons of 2, o-dibromoaniline , v.diich was studied as a 20 per 
cent solution in 001^, is fo-und to be predominantly due to 
external random fields ha-ving mean square value t-s'ri.ce at 
the site of ' proton compared to either of the two B 
protons. The random fields at the two 'b'^ protons are found 
to be highly correlated, vrhile those at '‘A^ and an.y of the 
'b^ protons being almost uncorrelated. Intermolecular dipolar 
interaction between the protons and the chlorine nuclei of 
the solvent appears to be the source for these random fields. 
Tho effect of intramolecular dipolar interactions between the 
various protons, ¥ri.thin the spin system, is small compared to 
that due to the intermolecular random fields. 


The relaxation of the AB spin system of the protons 
of 2-bromothiazole, studied as a neat liquid, is, on the other 
hand, found predominantly by tho intramolecular dipolar 
interaction between the two protons and by their scalar 
coupling with the nearby which inlfeum is relaxing fast 



173 


due to quadrupolar relaxation, iiie intramolecular dipolar 

interaction is significant in this molecule "because tlie 

proton-proton 'distance is small conpareato^in molecules like 

"befizene. External random fields do not seem to contribute 

significantly in this molecule. Double resonance experiments 

performed over a mde range of strengths of irradiation and 

at various frequencies of irradiation led to the determina- 

14 

tion of the magnitudns. and relative sign of profeon-H" 

scalar coupling constants. Ihc various parameters describing 

—12 

the proton relaxation are, = 7.0 x 10 sec,, •qs 9.7 

Hz and Jgjj ::r- 2,9 Hz (for = 0.61 millisec); and 

have the same sign. 

It may be noted that the presence of bromine does 
not modulate proton relaxation in both the above cases. The 
random fields arising die to th.. scale,r coupling of protons 
and bromine are negligible presumably due to a fast quadrupole 
relaxation of brominenuclei which possuss a largo quadrupole 
moment. 

Eurthor studies of this kind can be aimed at (i) 
studying the relative importance of dipole-dipole interaction 
and external random fields as a function of distance between 
the interacting spins, (ii) studies by double resonance in 
which scalar coupling. is only partly “ washed-out"^, (iii) relaxa- 
tion effects in double resonance as a function of external 
parameters like temperature, external magnetic field etc., 
for determining the relative inportance of various 
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interactions, in particular, that of the intramolecular and 
intermolecular dipole-dipole interactions, spin-rotation 
interaction, anisotropic chemical shift (if any)* Further 
insight into the relaxation processes can also be gained by 
attempting to interpret the different values of the mean 
square random fields and their mutual correlations in terms 
of specific models for the djoasmical behaviour of the 
li quid. 
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LISTING 


APPENDIX I 
OP COlilPUTSR PROGRAM! 
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rJTRODUGTIOIT 

‘The listing of a typical conputer program used in the 
density-matrix analysis of the double resonance spectra, is 
given in the follomng. The prograra is written in Fortran IV 
and run on IBM 7044. This program in its present form is 
applicable upto three spin systems, but can easily be extended 
to larger spin systems. Further, this program calculates doijble 

, I 

resonance spectra for high values of v^, for which Bloch 
approximation holds. For ^ Ioyj ' strengths of irradiation, the 
laboratory coordinate calculations have been performed by 
incorporating minor modifications to include the off-diagonal 
matrix element of IC . 

In the first part of this program the spin hamiltonian 
matrix is formed in the spin product basis. A matrix dia- 

gonalization routine is used to obtain the eigenvalues.”^ This 
routine also gives the transformation coefficients which are 
then used to transform various quantities from spin product 
basis to a,a'^.o. basis, in vAiich is diagonal. 

In the second part of the program, various relaxation 
matrix elements (R) are calculated by using a subroutine named 
RFDGIL. The relaxation matrix elements can be calculated for 
any degree of correlation. The simultaneous equations for 
matrix elements of X- are then formed and solved by using a 
determinant evaluation subroutine named BSTIR. In the last 
part the relative intensities of the double resonance trsmsi- 
tions including relaxation effects are calculated. 

+ This routine and a part of the remaining program are similar 
to those used by the magnetic resonance group at the Depart- 
ment of Ghemistry, Harvard University, Gambridge, Mass., 
U.S.A. during 1961-65. 
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SJCs Ti:'E008, P';r-ESf'30, MAME ANIL 

S T3JCE 

. ZdFTC yATl CDECK 

C TLi:. p.- ,-:GR.-a;' CAlCoLATEc DOUolE RESONANCE SPECTRA uSiNG EXTRE;iE 
V '•'m"' 'OL IR6 /^iS'JO.PT IC*' AMO -LOCH APPROX T ’AT I ON . 

C r.'SPlr- IS THE C?0 £r aF ThE SPIp SYSTE.'-. COKSI DErED, NASSYM IS THE 
C HHL^.EP JF ASSY'IETRIC STATES. CHART IS THE SPEED OF CHART PAPER 
C Pi_F: C.poS. GIVES DIoECT CmLI-ERATICM OF THE SPECTRA. 

C THE F ARE CoaPGNAMTS OF LATTICE FUi'iCTIONS IN THE IRREDuCIBIlE FORM 
C THESE riAY BE EXPRESSED LITHIN A PRO'^ORT I ORAL I TY CONSTANT 
C -.X IS THE NUMBER OF CO 'PCNANTS OF LATTICE OPERATORS. 

C AS STORES THE MATRIX ELE.'-.triTS CF THE RELAXATIONS HA[ s I lTON I AM I, ns THE 
C “^dI t product bASISs FED AS ImPuT DATA. THE CHEMICAL SHIFT C A-b 5 , 

c Ti E coupling Constant r» the amplitude of irradiation va=Vd and the 
C FrELsUEnCY 01 IRRADIATIO' (A = FRtQ.0F A RESONAmCE- FrEq* OF IRRAD. 

C Amu u-FrEq. oF B rESOI.An’CE - FrEq. oF IrrAD.) ArE INPUT PARAMETERS 
C Ul, C2» and C3 are the DEGrEE OF CORRELATION AND ArE TO BE FED AS 
C I! PUT PAFAi'SETERS. 

C N D R FOR BLOCH APPOCX TIAT I ON 
Mi'iDR THREE SPIN 1/2 

Di/iENSION ABC(3 ,S,15 ),F(15),AS(3,8»15) ,R00T( 10 > ,AA( 10»10) 

OIME.'.SION X(6»8),Y(d,8),E(8-.8),H(8) 

DIMENSION UdO.lO), V ( 8 ) , S I GMA { 10 , 10 ) jREPV EC ( 10 ) 
dimension PAAbudO.lO) sPADAb{S»8) ,RAA3C ( 8 . 8 ,8 ) 

COM', ON NOSPIFS >NQMAX»F 
CCr'.’-T/'! /ANIL/ABC.Cl »C2,C3 
com; on / KUMAR / ns YM 
STV,'0=SQRT(2.r ) 

READiOO? iCHART 
PRiNTlOlO.CHART 

A HF RD10u 9 ,NOSP I U, NASSYM 
NDi-; = 2^-i:-HOSPlN 
ISY.' =NDIi''i-NASSYM 
PR I NT 1006 

7 READ1C05 
PRIN.TIOOD 
REAu1009,FsqMAX 
DO 12 sNQ=l ,NGI'';AX 
oo i2C I=1,NDIM 
DO 120 J=l>NDIi'^ 

12c ASd »J.N0)=C.C 
D0121NQ=1,NQ iAX 
11^ READ1C30, I ,J>ASS 

IF( I oEQ.Oigo to 121 
AS( I ♦J,NQ)=ASS 
PPINTIOSI, r ,J ,:-.1Q,aSS 
GO Toil 9 
121 CONTINUE 

8 READ1007. (Fd ) »I=1»M0MAX) 

PRINT1007, CF{ I ) d = l*NaMAX) 

C 

C Xd.J) are THE matrix ElEMENTS OF THE SPIN HAt«ILTCNlAN IN THE SPIN 

C I RODUCT BAS^S ^ 

C FOR DIFFErEi.T SPIN SYSTEMS THE FOLLOWING CARDS HAVE TO BE CHANGED 

C THE PRESENT SET OF CARDS ARE FOR AN A82 TYPE OF SPIN SYSTEM. 

C THE ANTISYsMMETRIC STATES, IF AMY, SHOULD BE THE LAST PAIR OF STATES, 



i lOCTjAsrjRilVAjVb 

PR i'''Tl'S06 

'^TlO' .}. A,B»R»VA,Ve. 

C CO' PCTE X ATRIX 

^ DO 5i = 1,MDI *! 

’''05.1 = 1 
5 X(I,J)=0. 

X( 1.»C) = (A‘{'2o'*^B'*"R)/2o 
X ( 2 )2 ) = (-A+2„-;'-b-R ) / 2o 
X{3,3)=A/2, 

X{4»4)=-a/2, 

X(5,5)=(A-2c*B-R)/2. 

X ( 6 » 6 ) = { -A-2 « *B+R ) / 2 . 

X{7»7)=A/2„ 

X(8,8)=~A/2« 

X( l,3)=V5/STl'iO 
X{ l,2)=VA/2, 

X(2,4}=VL/ST',.fO 

X{3,5)=VB/STOO 

X(4,6)=VB/STW0 

X{3,^)=VA/2. 

X(4,5)=R/5.TU0 
X(2.3)=R/5TW0 
X( 5,6)=Va/2. 

X(7,fa}=VA/2. 

DO 3 1 = 1 ♦NDP^ 

DO 3 J= I, HO I'M 

3 x( j, n=x( I , j) 

c 

C '.kTRIX diagchalization rout I, he 

C ,-.RRAYS USED X(I»J)» E ( I , J ) ( E i GNVECTORS 3 » D I Y>ENS I OH-NDIM 

C FI, HU LARGEST OFF 01 AGONAL X 

C 

25j DO 2b3 I = l,HDI.'l 
00 253 J = l,\!Dir! 

IF(I-J) 251,252,251 

251 E(J,J)=0. 

GO TO 253 

252 r{l,J)=l,0 

253 COHTINUE 
2c^ LA=1 

UP = 2 

UX=ABS{X(2,1) ) 

DO 2u2 J=3,MDIM 
JA=U-1 

DO ^u2 1 = 1, JA 
0=1X“ARS(X{ I ,J) ) 

IFfO) 2''1,202,202 

201 r2X = ABS<X( I ,J) ) 

LA=I 

LB = J 

202 CONTINUE 

IF (BX-O.OOOOl ) 320,320,300 
C COr’PUTE ROTATION 
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3‘ • i:x = A ( LA yLc ) -"-X ( LA > LB ) 

D=X(LA,Lt)-X(L5»LB) 

.-^=SL'5i (mFS( (';:+D)/(2«c-»-q) ) ) 

IF (<- oTOT-A) 302,3o2>3Q1 
3-1 3 = -'. ''-A 

A=S33T ( lo 0-6*3) 

CO TC 303 

3^2 l=-S':RT { 1.0-A*A) 

3».3 IF(L'/X(LA 5LB) } 3CA » 3 lO , 3 1>J 
3' 4 P.= C.0-F 

C GRTHOGC^ ''AL ROTATIOr: OF '^lATRIX 

Sir. r>o 313 J=1,MDI.“ 

IF(LA-J) 311,313,311 
31] IP{Lr:-J) 312,313,312 
3ii X ' lA , J ) =A*X ( J ,LA )-B*X ( J , LB ) 

X i Lb » J ) =B*X ( J , LA ) +A*X ( J » LB ) 

313 CONTINUE 

no 314 J = l,MDl!''f 
X( J»LA)=X(LA,J) 

31-4 X{ J,LB)=X(L3,J) 

r)=X(LA,LA}+X(LB,LB) 

X{LA,LA)=A*A*X(LA,LA)+B*B*X{L3,Ld)-2 

X(L&»LR )=D-X{1 A, LA) 

X(LA,LB}=0.0 

X{LF,LA}=0,C 

C Ef'^’ER ROTATION IN E i'ATRIX 

no 315 J=1,NDIN 


*A*B*X ( LAsLB) 


5 X = E( J,LA) 

E ( J 5 LA ) = A*DX-B*E ( J , LB ) 

31 j E( J»LB)=B*BX+A*E{ J,LB) 

GO 0 200 

c TEST DIAGONALIZED MATRIX 

32' uX-.'.0D'04 
JA = '-''r 1 1-1 
no 223 1 = 1 , JA 
J b — I + 1 

no 323 J=JB»NDIM 
IF(X{I,J)) 321 » 323,321 
32 i R=AbG(X(I>J)/(X(I,I)-X(J»J) >) 

IF(K-BX) 323 , 323,322 
32 ? BX=R 
LA=I 
LB = J 

323 CONTINUE 

IF IBX-O . 00005 } 45 o» 450,300 _ 

C GALC TRANSITION INTENSITIES AND FREQUENCIE. 

DO 451 I = 1 »NDIM 
ri( I )=X( I ‘I ) 

00 451 J=i,NniNi 
451 X(I?J)= 0.0 


C 

c 

c 


TE FOLLOWING ArE MATRIX ELLf'.EXtTS OF H + J IN SPIN BASIS. 

-ESE cards lave to BE CHANGED FOR DIFFERENT SPIN SYSTEMS. 
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v( i,2)=l.r 
X ( 1 ,3 ) =STV;0 
>:i2,4)=ST'-.'0 

X( 3,5)=STV'0 
a{A-, 6)=ST .0 
A ( 5 » 6 ) = i c 
X{7»o)=l. 

432 DO ^53 I = 

DO 433 0=1, Nr I"! 

YU »J)=0«0 
DO “^‘->3 LA=1»NDIM 
Dr 453 l8=3»NDIM 

4'=^'3 Y ( i , J }=Y( I , i)+ E(LA,I )*X(LA9LB)-''-F(L[i,J) 

’Rli'JT 1001 
DO 456 I = 1,^01.''’ 

PRIf'T 1G62»I.H(I) 

DO 456 J = 1.MDF'; 

R=Hl I)-H( J) 

Y60= '{ I , J)*Y( I »J) 

PLOT = F '5-C‘^4RT 

456 PRINT lo(’3,J,R»Y5Q»Y( I » J) , hi I, J), PLOT 

C CALCULATIC.i OF REDFIELD ELE.-'EnTS 

C TPA .SFOR;' RELAXATIO"'' MATRIX 

D013 jfiO=l»NQMAX 
DO 130 I=1,NDI.':1 
DO 13o J = l,NDr''' 

APC( I ,J»NQ)=0.0 
nOl30K=l •NDIK 
DO 13n L=l,NDIi'l 

l3o ABC ( I » J »NQ ) =AbC ( I , J , NO ) +E ( < » 1) *AS ( K. » L »NQ 3 ^!-F { L , J ) 

C F:R,M SIGi'lA .'lATRlX At'^D TR Al'iSFOR^'l IT 
00515 1=1 »MDIM 
C0515j = l ,NDIR- 

5a3 X( I ,J)=0.0 

C THE FOLLOWING GIVE -F ( Z ) =-S'JO ( U Z ) ) IN THE SPIN PRODUCT BASIS. 

C THESE are thus MATRIX EbEr'ENTS OF SIGMA ZERO IN SPIN PRODUCT BASIS. 

C THEoE cards have Tq BE CHANGED FOR DIFFERENT SPIN SYSTEMS. 

C 

Xf 1,1)=-1.5 
Xt2»2)=-0.5 
X{3,3)=-0.5 
X(4 j4)= >.',5 
X{5,5)= 0.5 
X(&,6)= 1.5 
X(7,7)=-0.5 
X(8»8)= 0.5 

c TFA'ISFORM SIGMA t ATRIX 
PRINTIOZO 
D05i7i=l,NDlM 
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DP016J = 1 Jivrir’ 

I U ^ ‘I A ( I * J ) — 0 » 

'i05i6L = l 

005 j &K = 1 ! 

516 oIG'.m( I , J)=SIG.-!A( I , J ) +E { < , I) *X { K » L ) { L , J ) 

517 PRI'iTlUjas (SIGr''A( I ,J) » J= 1 , .VO I P' ) 

FPliiTl006 

r RFAo correlator 

6 RE^'O 1 )02 ,NCCRLN 
D0900 lJr=l,NCORLV 
C Cl, C2 AMD C3 ARE CORREi-AT I ONS . 

9 READi007,Cl,C2,C3 
45P PPira 1006 

PRINT1.C.2A,C1,C2 ,C3 
PRli‘STlC06 

C CALCULATICM of R(IIJJ) 

D0 51C 1 = 1 »NSYN', 

D0509J=I ,t\!SYM 

CALL RFOCAL(I»I,J,J»RFO) 

Raa3B( J,I )=RFD 
509 KAAOE( I ,J)=RFD 
COMTINUE 
D0512J=1 ,MSY.’''! 

D0512K = J,MSYi7 
IF(JoEQ.R) GO TO 512 
005 11 1=1 »MSYM 
CALL RFDCAL( I , I , J,K,RFD) 

T AAt'.C( T ,J,K) = RFD 

511 n/ AFC( I ,K,J)=RFD 

512 CONTINUE 
0051Ai=l ,NDIM 
D0515J=I ,140171 

IF ( I .EQ® J) GO TO 513 
CALL RFDCALI I , J, I ,J,RFO} 

RABAB( I ,J)=RFD 
RABAGI J, I )=RFD 

513 CONTINUE 
I I = T + 1 

51A PRlNTlvlS, C (RAdAB( I , J ) , I , J , I , J ) , J= I I , NDI M ) 
C SOLUTOION OF SI-IULTANEOUG EQUATIONS 
DOj20I = 1 ,NlviM 
/( I ) = 0. 

D0520J = 1 ,MSY!V. 

D0520K=J,NSYM 
IFfJ.EQ.iOGO TO 520 

v( I j=V{I )+2,0*SIGMA{ JfKI-s-RAABCI I,J»K) 

52r> CONTINUE 
PRii-iTl )18 
D052Al=l ,NSYM 

5 24 PRIflTlOlO , (RAABBI I »J) ,J=1»NSYM) ,VU ) 
PRINT1006 

CALL DETER(RAABB,DET ) 
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-F If.;T10ll ,DET 

L‘-‘,TCK PC'? LT'EAR DEPE'JCE''CE 
D0530K = 1 j.'iSYH 
D0 5i5l = l sf'iSY-'; 

D05Z5J = 1 .'YSY.M 
-cpv7C( I )=V( I } 

523 U{ I »J)=RAABE{ I >J) 
n'05<i6J=l ,f,SY-( 

REPVrC{!<.)=0»0 
CALL 0ETFF(J,DET) 

PRr''Tlr'2x,K,DET 
Ci FCK = ABS{DrT } 

IF ( C--ECK-O.C001 } 530 » 54:., 54; 

530 COMTIIsiUF 

54r IF (ChiECK.LT.C.OC-Ol) GO TO 900 
Of 5b..iA=l *FSYt'; 
no5TlJ = l ,fiSYY. 

D0571 1 = 1 ,f-iSY-l 
AA(I*J)=U{I,J) 

IFfJ-K) 571,570,571 
57o AA( I ,L)=REPVrC( I ) 

3 71 CC,\TIMUE 

CALL DFTER(AA,DEM) 
rvOOTCK.) = PEFi/DET 
PR I Ml TOS 

58r PRlMlnl2,f<,DEi-i,K,R00T(K) 

C CHFrc THF ROOTS 
PRinTinC6 
’Rr'TlCl6 
005E4K=1 ,M5YM 
EG=n.o 

D0533 J=1 ,NSYM 
563 EO = £:.''+ROOT( J)4U(i<,J) 

PSh PKli."1017,K:,EQ»REPVEC{i<.) 

C CALCULATION OF SIGNAL INTENCITIES 
rRiNTlnof 
PRlOTln26 
0061;. ! = 1 ,NDIM 
PRlMTlO-^Z, I ,H( I ) 

OC610J=1 ,NDIM 

JF{ T oEQ» J) GO TO 610 

SI G=Y{ I , J ) *Y{ I » J )*( SIGMA ( J , J ) -SIGMA (1,1 >+RO0T { J) -ROOT ( I ) ) 
SIGI NT=STG/Ra5AB (1,0) 

FRE9=H( I )-H( J) 
plot=frfq-»chart 

PRINT10C3,J,FREQ,SIG,SIGINT,PLOT 

610 continue 

9-^n CONTINUE 
PRiNTlOOf 
PRlNTlnOb 
GO TO 1 
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iC r( s-.-^vT ( A=,^lG.5»4h C=,FlC!o5»4H P=,FlC.5»5H VA=,F7»3,5h VB = 

l»F7e3//^ 

iuol rCF AT [2ri 1 ,6X ,4Ht ( I ) ,6X »2ri J , BX , 4HFREQ , QX » 3 H I WT , 8X , IQHTR AiXS (1+) 
1,4X511HTRA!MS C0EFF»* CHART^/3 

i003 FOR -'.at (I2,F12o 5) 
i."-'-' format (2':x>I2»5f14.5) 

IpfA FORi’iAT (2F1S.& ) 

IC'b FCRF'AK-';- 
iQ^'o F0R7AT(/) 

3cG 7 F0';''-.AT(BF10<,5 ) 
lOGd FOR'-’AT { 10F12 .5 } 

FOG-;AT(5i5' 

miu FCR-.AT( * ■'iO OF CMS. PER Cop.S. =-)«■ , F 10 c 5 / ) 
loll FOR-iAT DET = -» ,El2 . 5/ ) 

1U12 For..AT{5H DEFi{»I 1,4H) = ,El2.5»lOH R00T(,I1,4H) = jElEoE) 

lulP FORFAKACFIS.B^* = R(*»All,iH*) ) 

1016 FoR'^AKEPH CAlCULAThO VALUE OF i_QUAT ION » 5X » 15H SUPPLIED VALUE) 

1017 FOR ‘AT ( 5X,1H( » I 2 ,3H ) = , FIG .6 » 1 5X > FIO » 6 ) 

1018 FGR:,mT(4uX»*A< I » J )*»50X ,-x-V( I )*/) 
lol^ FOPFiAT (8F12.6,6X »F12.5 ) 

1020 FORl'iAK //35X» ■’-SIGMA ZERO MATRIX-'*-/) 

1021 FOP. AT{-» EOUATION REPLACED = * » I 2 DET = -”- , E 12 . 5 / ) 

x 022 FOPy.AT( 1CX»* MATRIX ELEMENTS OF !( + )-;;-) 

102' F0R"iAT<2l5»*=*»F10.5 ) 

1024 FOR^’AK'!*' C1 = -5^»F10»5>-!''- C2=-k-»F10o 5»^ C3=-i*'»FlO*5) 

1026 F0R-'AT(-«- I-it »6X 5 *E ( I ) * > 8X J* , 8X » ->FREQ-“- s 1 IX »■**■ I NT-“- > 6X > * I NT /R I J I J* » 

18X»*CHARTr-/ ) 
i03,; F0RMAT{2 t5,F13.3) 

1031 FORmAT{5< AS(*,I1,-:;-,-k-,I1»-"-,*,I2,4}='S*-,F13„8) 

END 

SIE^-'TC SUEl i'iODECK 

SUBROUTINE RFDC AL ( I » J » K » L » RFD ) 

C '■hIS subroutine calculates THE RELAX AT T O!" tMATRiX ELEMENTS Ro 
C I ''.RTIAL CO^.r ELATION WITH EITHER EXTiB'ISS OF COfPLETF OR NO COREEL AT lOI 
C IS possible, the INPUT IS FROM .')AP1 PROGRAM AI'D THE RESULT IS 

c storfd in RFD. 

DI -lENSlCM A( 8 >3, 15 ) ^F ( 15 ) 

COM ;N N0SPI11»:''!C’'iAX,F 
COf. GN/ANIL/ A > Cl»C2,C3 
f4 I ' =2^c«-KOSPlN 
. 00P=NQ,-1AX/N0SPIN 
RFD=o<,n 

DO 230 NQ=1,NQMAX 
T+NGOP 

I :•=. .3+noop 

LC = !’.0-N00P 
JQ=LO-NOOP 
RQ=OoO 

RQ=.^0+A{ I »K»NQ)-»-A{ J»L,iMQ)-!fF(;'iQ) 

I F ( f-Q .LE . NOOP ) RO = RQ-^-Cl-» SQR T ( F ( NO ) *F ( HQ ) ) *A ( I » K » NO ) *A ( J » L , MQ ) + 

1C2 "-SORT I F ( ^'0) -^-F ( I Q ) ) *A ( I » K » NO ) *A ( J » L , 10 ) 

I F ( NO . GT . KUOP . AND . L.5 o LE . NuOP ) RO= RQ+C 3-“-SQR T ( F ( NO ) *F ( MQ ) ) * A { I » K ,.NO ) 



n n n 
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2 x: 
211 
2 '’ 2 


21 


2 j.<f 

215 


216 

21 I 

220 

221 

222 

223 

224 


1^* ( J »L 2 )+Cl->3-.RT(F ( LO) "-F ( NC ) )-;fA ( I )*A ( J ,L»LG ) 

IF ( J. oGT, ' 0:.=^) R j = R-:v+C2-x-a ( I ,k , 0 )*A ( J,L, JC )"^SCRT (F(f';Q)"'«-F ( JQ) )+C3rr 
ibc-'^T ( F( ';g)-'--f { lG ) )*A( I {J ,L»L0) 

IF(J-L)22r:,210>220 
IF{ ''0-F'?P)212»212»211 
IF(L2 --v:'0P) 214,214,216 
G0213 :'!=! 

'<C = ' 0— A ( f » K 5 ' ■ ; ) A { , I » I , ■ 0 ) *F ( 'Or) ) /2 , 0 

.■>C = ;Q-C1-"-a {.' »< 5 'C ) 5 1 ,i'''G)*5CRT ( F ( WQ )^-F(r‘a ) ) /2oO-C2-''-SOf5T { F ( NQ) * 

1 F{ !u) )*A(,’-;,K»f'-..-)-A{.;»Ii.IG)/2.0 
GO TO 22n 
D0215 m=1 ,NDIW 

Rr = r 0-A < I ‘ » K »M0 ) *A { ■ 1 » I 9 . ’0 ) *F I • .0 ) / 2 . 0 

.M =kO-C3*A(tM9K»KQ)--A (;;» I 9 .0 ) -“'SC R T { F ( MQ ) ’“-F ( f '0 ) ) / 2 . 0-C l*SGRT { F{ I'iO) * 
IF(LO) )*A(,M»K9Na)"fA(.'.9 1 9LQ)/2,0 
GO TO 220 
002i7M=l ,NDIM 

RQ = U 3-A ( V 9 K 9 r;Q ) ■» A ( 'm 9 1 9 A'O ) -“'■F ( GO ) / 2 , 0 

^ 0 = ;- 0-C2 *A ( V] , K 9 AG ) * A ( , I , JQ ) ---SORT ( F ( MQ ) *F ( JQ } ) / 2 . 0-C 3*S0R T { F ( NO ) * 
IF ( LO ) ) *A { ;Vi 9 X 9 NO ) -:f A ( f'. 9 1 9 LO ) / 2 » 0 
IF( I-X)230,223 ,23n 
IF ( f:0-!'.OOP) 223 ,223 ,2 22 
IF(LQ-POOP)225,225,227 
D02 24 -1 = 1 ,NDIM 

r0 = A 0- A ( M 9 i 9 , GO ) ( M 9 J 9 0 ) F ( NO ) /2 ,0 

RQ = RC-CHf-A( l9L9iJG)-:fA( i,J, ;G)*5GRT(F(NO)-'^F<:;Q)}/2.0-C2*SORTCF(NQ)^«- 

li- ( IC) )*A(;'9L9*'0)*A(M9j9lQ) / 2 o 0 


GO TO 23 0 
225 D022fe[- = 1 ,,MDIG 

RQ = r.:;_A (,-1,2 9 NO ) ( . ! 9 J 9 •’ 'C ) F ( ,■ 0 ) /2 „ 0 

2 26 ro = ,iO-C3*A{'.9L9iVn)^:-AC'i9J9: 0 ) -'^SGPT ( F ( MO ) *F ( GQ ) )/2.0-Cl*SQRT(F{NQ)* 
IF ( LC j ) -' 1 - A ( ; ’ 9 L 9 riQ ) ^:-A ( .'■•'i 9 J , LQ ) / 2 o G 
CO TO 230 
227 00228 m = 1 9NDIi'i 

R0=R 0-A ( i-t 9 L 9,00 } -"-A { i‘t , J , GO ) *F ( !- 0 ) / 2 „ 0 
22 b :nO = .^C-C2*A(i'UL 9N0)*A(.'i, J, JO) --'SORT ( F (MQ)*F( JQ) ) /2,0-C3*SGRT { F( NO) * 
IF ( LG 3 ) *A { M 9 L 9 NO ) *A ( I'i 9 J 9 LQ ) /2 o 0 
23o RRr)=RFD+RQ 
RFT'JRM 


FMO 

6I3FTC SU32 MOOECtC 

SUBROUTINE DETER (B 9 DET) 

This subroutine calculates THE VALUES OF A DETERMINANT BY PIVOTAL 
HETHoD, NO IS THt ORDEr oF THE DETERMINANT TO BE CALCULATED AND 
DET IS ITS CALCULATED VALUE. 


COY,- Crj/KUr-AR/NO 

DiHEMSiOf 8(10910) ,A(10,l0) 9pR0U( 10) 

l^u^'l=^:o 

DFT=1. 

D03I=1,N0 



f>C3J = l ,M0 

3 a( I ,j)=a( I , J) 

4 j06 

IF C A ( 1 ,y,)-0o )9 »6 ,9 

6 CONTINUE 
DET=0o 
GO TO 18 

9 PwOr=Ml,M) 

DO 7 I=i,\UH 
f^ROW( I ) = ^ ( I »K} 

7 A<l,I)=A(l,I)/PlVOT 
DO 13 i = 2,NUi'''! 

D013 J=1,NUK 

13 AU »J)=An ,J>-PROVJ{ I)*A{ 1»J> 
D020I=2,NUM 
D020J = 1 »NUr''' 

IF( ,20»1C 

10 A( I-1,J-1) =A{ I »J ) 

GO TO 20 

11 A{ I-lsJ)=A{ I ♦J) 

20 CONTINUE 

NUM=NUM-1 

DET = PIVO^^'<-OET-'!- ( (-1. ) •"■*( >1+1) ) 

IF{NU:T-2 ) 18,16 

16 DET=(AC1,1)*A(2,2)-A(1»2)*A(2,13 )*DET 

18 return 

19 END 
SFNTFY 
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and Solvent Effects on Coupling Constants. 


Heprintcd fron llolecular Physics, 1967 
Vol. 12, lTo.6, 593-596. 



AII-2 


Reprinted from: 

MOLECULAR RHLSICS, 1S67, VCLi. 12, No. 6, 593-596 


II.M.R. SPECTRA OE 1 , 5-I)LPLUORO-4, 6-DINITROBENZEEE AND 
SOLVENT EEEECTS OIT COUPLING CONSTANTS 

oy AITIL KUMAR 

Department of Physics, Indian Institute of Technology, 

Kanpur, India 

(Recieved 13 Eehruary 1967) 

1. INTRODUCTION 

The high resolution nuclear magnetic resonance 

spectra and the solvent and concentration dependence of the 

spin-spin coupling constant of 1,3-difluoro 4y6-dinitro'ben- 

zene are reported here. This molecule is a four spin system 

of the type AjBX 2 notation of [1]) with nuclei 

labelled as shown in figure 1. 

(B) 

H 

(X)E-^^^XE(X) 

.1 i : 

H 

(A) 

Pigure 1. l,3-difluoro~4, 6-dinit rohenzene. 
Analytical expressions can he derived for all the allowed 

r , 1 

transitions and their relative intensities L2J. The H 
and ^^P resonance spectra of the molecule, recorded on a 
Varian HE- 100 spectrometer operating at 100 Mc/s and 94*1 
Mc/s respectively are shown in figure 2. The observed 




- 4 . 1 1 — 5 - 


5’ig. All- 


2.i!'luorine , X proton, ond B, spectra of 
1, 3-difluoro-4, G-dinltrobeiizone. file lines 
in A resonance are broader compared to B 
resonance, probably due to a smaJLl loupling 
of the A proton v/ith nijrogen, in the nearby 
NO 2 group, which is relaxing rapidly due to 
quadrupole interaction, fhe separation bet- 
ween the outer lanes of the tripclots in A 
and B resonances is 2'Jp-p rud respectiTely. 

In X resonance the separation betvmon the 
outer lines is ( 1 + j Jppj ) that 

between the inner lines is (|jp.i;,, j - 



4 
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spectra, liov;ever, do not shov.* any further structure, 
indicating that the para H-H coupling is less than 0.3 c.p.s, 
(^half the linevh.dth) , Thus if is assumed to be zero 

for this system the observed proton and fluorine spectra give 
in a straightfor\7ard marmer the values of the meta and 
ortho H-? coupling constants. 

2, SOLVENT EEEEGTS 

1, 3“ difluoro-4, 6-dinitrobenzene is a solid at 
room temperature and can be studied in high resolution 
through the use of a solvent. The spectra were studied in 
various solvents to see the effect, if any, on the spin 
coupling constants due to the solvent. The results are 
given in the table. In two of the solvents deutero- 
chloroform and dioxane the spectra were studied at several 
different concentrations. The variation of the coupling 
constants mth concentration, in the range of concentrations 
studied, is smaller than the variation in different solvents* 
The coupling constants tend to increase vdth percentage 
concentration in d-chloroform. This trend is reversed in 
dioxane and further the values are smaller in d-chloroform 
than in dioxane. The coupling constants in the pure 
compound therefore are expected to lie between the values 
in these two solvents. 

The solvent and concentration dependence of 
coupling constants, except where they are due to configura- 
tional modifications, indicate changes in intramolecular 
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interactions, miie various factors considered responsible 
for these changes in and coupling 

constants, reported hy several workers [3-8], include, the 
reaction field effects of the medium, molecular association, 
complex formation and the shape of the solvent and solute 
molecules, Ihe reaction field [9] is an electric field 
produced at the solute molecule due to a polarization of 
the medium hy the polar solute. She effect of this is 
indicated hy similar changes in chemical shifts and coupling 
constants with dielectric constant of the solvent [3-8], 

The present observed variation (table) of the sum 
of the two coupling constants 

cent concentration in the first four solvents is approxi- 
mately linear mth dielectric constant. This variation is 
therefore likely to be due to the reaction field effects of 
the medium [3-8], In the other three solvents namely 
acetonitrile, tetrahydrofuran and dioxane, the coupling 
constants bear no such relation to the dielectric constants. 
These three later solvents are kaown to be high donor centres 
[10] and will probably form donor-acceptor complexes by 


This sum, which can also be obtained directly from the 
fluorine spectra, shows a larger variation between solvents 

O HI 

than either or Jttt^ separately. The similar 


HF 

variations of J 


o HF 

and Jttti with solvents indicate that 


HF HF 

these two coupling constants are of the same sign [6]. 
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hydrogen bonding between the ring protons of the solute 
molecule and the donor centres of the solvent molecules 
[4-, 7]. Ihe complex formation may affect the coupling 
constants in a v/ay which is different from the reaction 
field. 

Both the reaction field and the complex formation 
are also kno\m to affectthe chemical shifts [8,11] and 
therefore it is reasonable to expect some correlation 
between the variation of the coupling constants and that 
of the chemical shifts. An approximate correlation is 
obtained in the present case. The two coupling constants 
tend to increase v/ith increasing chemical shifts, and 

Vg and decreasing values of ( in agreement with 

previous observations [3-6]. The chemical shifts in 
benzene are, however, very different from those in other 
solvents, probably because benzene is an aromatic compound. 

These observations suggest that the effect of the 
solvent on the coupling constants should be taken into 
account, particularly in solid samples where solvents are 
unavoidable. The apparent concentration independence of 
coupling constants within experimental errors could, in 
several cases, be due to the limited range of concentrations 
v/hich could be studied (limited by solubility, signal-to- 
noise ratio, etc.) rather than the absence of any solvent 
effects. Therefore, the solvent dependence of coupling 
constants should be determined by studying the concentration 
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dependence in ai least two solvents of different natures, 
and tile pure compound values obtained by extrapolation, 
oust as is tbe practice for chemical shifts. Conclusions 
arrived at otherwise could be erroneous [12], She study 
of solvent dependence of the coupling constants should also 
throw some light on the various intramolecular interactions 
responsible for the coupling between two spins in a ' 
molecule. 

The author wishes to aclmowledge many invaluable 
discussions with Dr, B.D, Hageswara Rao and the keen 
interest taken by Professor P. Yenkateswarlu in this work. 
Thanks are also due to Dr. Ramesh Chand for sending the 
sample from the U.S.A, 
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